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> '■ ABSTRACT 

' Using Far Ultraviolet Spectroscopic Explorer (FUSE) and Hubble Space Telescope 

(HST) observations of the QSO PG 1259+593, we detect D I Lyman-series absorption 
^ ■ in high velocity cloud Complex C, a low-metallicity gas cloud falling onto the Milky Way. 

. This is the first detection of atomic deuterium in the local universe in a location other 

than the nearby regions of the Galactic disk. We construct a velocity model for the sight 
line based on the numerous O I absorption lines detected in the ultraviolet spectra. We 
I identify 8 absorption-line components, two of which are associated with the high velocity 

■ gas m Complex C at « -128 and ^ -112 kms"^. A new Westerbork Synthesis Radio 

Telescope (WSRT) interferometer map of the H I 21 cm emission toward PG 1259+593 
indicates that the sight line passes through a compact concentration of neutral gas in 
Complex C. We use the WSRT data together with single-dish data from the Effelsberg 
. 100-meter radio telescope to estimate the H I column density of the high velocity gas 

^ ' and to constrain the velocity extents of the H I Lyman-series absorption components 
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observed by FUSE. We find N(H I) = (9.0±1.0) x lO^^ cm-^, N(D I) = (2.0±0.6) x 10^^ 
and N(0 I) = (7.2 ± 2.1) x 10^^ cm"^ for the Complex C gas (68% confidence 
intervals). The corresponding light-element abundance ratios are D/H = (2.2 ib 0.7) x 
10"^ 0/H = (8.0±2.5)xl0-^ and D/0 = 0.28±0.12. The metallicity of Complex C gas 
toward PG 1259+593 is approximately 1/6 solar, as inferred from the oxygen abundance 
[0/H] = — 0.79ibQ;Jg. While we cannot rule out a value of D/H similar to that found 
for the local ISM (i.e., D/H ~ 1.5 x 10~^), we can confidently exclude values as low as 
those determined recently for extended sight lines in the Galactic disk (D/H < 1 x 10~^). 
Combined with the sub-solar metallicity estimate and the low nitrogen abundance, this 
conclusion lends support to the hypothesis that Complex C is located outside the Milky 
Way, rather than inside in material recirculated between the Galactic disk and halo. The 
value of D /H for Complex C is consistent with the primordial abundance of deuterium 
inferred from recent Wilkinson Microwave Anisotropy Probe (WMAP) observations of 
the cosmic microwave background and simple chemical evolution models that predict 
the amount of deuterium astration as a function of metallicity. 

Subject headings: cosmology: observations - Galaxy: evolution - ISM: abundances - 
ISM: clouds - nucleosynthesis - quasars: individual (PG 1259-1-593) 

1. Introduction 

Observations of the abundance of deuterium relative to hydrogen (D/H) in different environ- 
ments provide insight into the evolution of the light elements in the universe. With the excellent 
concordance in the estimates of the cosmic baryon density from measurements of D/H in low- 
metallicity quasar absorption-line systems and measurements of the cosmic microwave background, 
the cosmic baryon density and the primordial value of D/H are now tightly constrained (Buries, 
Nolett, & Turner 2001; O'Meara et al. 2001; Spergel et al. 2003). Thus, it should be possible to 
test chemical evolution models by examining the progression of D/H with time (see Lemoine et al. 
1999 and Olive, Steigman, & Walker 2000 for recent discussions). The abundance of deuterium is 
expected to decrease with time since there are no known sources of deuterium capable of increasing 
the cosmic abundance significantly (Epstein, Lattimer, & Schramm 1976; see also Prodanovic &: 
Fields 2003). Many chemical evolution models predict moderate levels of deuterium destruction by 
stellar nucleosynthesis, typically less than a factor of 3-5 (Clayton 1985; Edmunds 1994; Steigman 
&; Tosi 1995; Tosi et al. 1998). More recent models suggest that slightly lower levels of astration 
are also possible (e.g., Chiappini, Renda, & Matteucci 2002). 

Unfortunately, measurements of D/H are particularly difficult, and there have been relatively 

few opportunities to directly measure the detailed changes in the abundance of deuterium as a 
function of metallicity or time. A key piece of information missing in discussions of the evolution of 
the light element abundances with time is the behavior of D/H in environments with metallicities 
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between those of the high-redshift systems (typicahy Z < O.OIZq) and those of gas in the local 
neighborhood of the Sun (typically Z ~ Z©). 

There are several reasons why it it is important to determine the D/H ratio in a wide variety 
of galactic and cxtragalactic environments. First, there have been few high-precision estimates 
of D/H at moderate to high redshifts (z ~ 2 — 4), where the amount of stellar processing of 
deuterium is presumably low, as evidenced by low metallicity (e.g., O'Meara et al. 2001; Pettini 
Sz Bowen 2001 and references therein). The measurements that have been made appear to yield 
conflicting values for the primordial abundance of deuterium, with the observed values depending 
on the type of system observed (e.g., Lyman-limit or damped Lya systems - see Pettini & Bowen 
2001). Second, estimates of D/H in many locations yield a more global perspective of the chemical 
history of gas at different epochs than is possible from a few isolated measurements. Chemical 
evolution models seeking to describe the general evolution of the light-element abundances need 
a large sample of measurements to avoid systematic problems encountered by relying upon data 
for only a few types of environments. Third, although measurements of deuterium in nearby gas 
clouds imply a relatively constant value of D/H within the local interstellar medium (ISM; Moos et 
al. 2002 and references therein), substantial variations may exist in D/H and D/0 over distances of 
only a few hundred parsccs (Jenkins et al. 1999; Sonneborn et al. 2000; Hoopcs et al. 2003; Hebrard 
&; Moos 2003). If a sufficiently large number of high-precision D/H and D/0 measurements can be 
made in a diverse set of nearby environments, it may be possible to understand the exact causes 
of this variability and the degree to which galactic chemical evolution and accretion of intragroup 
gas clouds influence the scatter in the observed ratios, both locally and at high redshift. This goal 
is a major science driver for the Far Ultraviolet Spectroscopic Explorer (FUSE) mission (Moos et 
al. 2000). 

To bridge the gap in D/H between low and high metallicity environments, we have obtained an 
extensive set of FUSE, Hubble Space Telescope (HST), and interferometric H I 21cm observations 
of the quasar PG 1259-1-593 behind high velocity cloud (HVC) Complex C. The HVC is located at 
least 3.5 kpc from the Galactic plane (Wakker 2001), well beyond all Milky Way clouds with current 
D/H determinations. Unlike previous investigations of D/H and D/0 in either high-redshift clouds 
or the local ISM, we know which gaseous system is responsible for the high-velocity D I Lyman- 
series absorption observed toward PG 1259-1-593. A global description of the gas in Complex C is 
available from both emission and absorption-line measurements. The neutral gas in Complex C 
has been mapped extensively in H I 21 cm emission (see Wakker 2001 and references therein) and 
low-ionization absorption (e.g., Wakker et al. 1999; Richter et al. 2001b; Collins, ShuU, & Giroux 
2003). The ionized gas in Complex C has been investigated in absorption by Sembach et al. (2003) 
and Fox et al. (2003), and in emission by Tufte, Reynolds, &; Haffner (1998) and Wakker et al. 
(1999). Complex C is an excellent site to determine D/H for comparisons with the high-redshift 
values because it is chemically young (Richter et al. 2001b; Collins et al. 2003; Tripp et al. 2003) 
and has a metallicity (10-25% solar) lower than that of the general ISM of the Milky Way and 
higher than that of intergalactic clouds at high redshifts. 
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In this paper we describe these new measurements and the resulting D/H and D/0 ratios 
in Complex C. In §2 wc describe the FUSE and HST Space Telescope Imaging Spectrograph 
(HST/STIS) absorption-Iinc observations and the H I 21cm interferometer observations. Section 3 
contains a short summary of the properties of the PG 1259+593 sight line. In §4 we outline the 
methods and general assumptions used to determine the column densities of H I, D I, and O I in 
Complex C. We determine the H I column density from interferometric H I 21 cm emission data and 
use the FUSE and HST/STIS ultraviolet absorption-line data to determine the O I and D I column 
densities. Sections 5, 6, and 7 contain descriptions of these determinations for oxygen, hydrogen, 
and deuterium, respectively, and provide estimates of the various errors. The column densities and 
error ranges are summarized and discussed in §8, and comments on future progress appear in §9. 
We summarize the results of the study in §10. 



2. Observations 

2.1. Far Ultraviolet Spectroscopic Explorer Observations 

We observed PG 1259+593 on nine occasions with FUSE between 2000 February 25 and 
2001 March 28 for a total (orbital day + night) integration time of ~ 600 kiloscconds (ks) spread 
over 236 individual exposures. Approximately 350 ks of data were obtained during orbital night. 
An observation log for the nine visits is provided in Table 1. Data were obtained through the large 
(LWRS; 30" x 30") apertures in all four FUSE channels (LiFl, LiF2, SiCl, SiC2) with varying de- 
grees of success. PG 1259+593 was always well-centered in the LWRS apcrtmc of the LiFl channel 
used for guiding, but thermal effects caused the light of PG 1259+593 to drift around inside (and 
at times even partially outside) the remaining three apertures. Detector high voltages were also 
at reduced levels during some of the exposures as a result of operational difficulties (see Table 1 
notes). We obtained all of the data in time-tagged photon-address mode to allow for data screen- 
ing and time-dependent corrections necessary to fully calibrate the data. The raw FUSE data for 
PG 1259+593 can be found in the Multi-Mission Archive at the Space Telescope Science Institute 
under the observation identifications P1080101-P1080109. 

We processed the data with a customized version of the FUSE pipeline software (CALFUSE 
v2.2.2), which is publicly available from the FUSE Project at the Johns Hopkins University.^^ The 
details of this processing follow the same general principles outlined by Sembach et al. (2001b). 
However, since PG 1259+593 is a faint object {Fx ~ 2 x 10~^^ erg cm~^ s~^ A~^ between 920 
and 1200 A), we modified these procedures in the following manner to process the data. First, we 
screened the raw photon lists in every exposure for Earth limb avoidance, South Atlantic Anomaly 
passage, pulse height distribution constraints, and particle event bursts (see Sahnow et al. 2000). We 
chose the pulse height restrictions after inspection of the pulse height profiles to reduce background 



^^See http: / /fuse .pha. jhu. edu. 
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events while minimizing the number of source events discarded; pulse height data numbers from 4 
to 24 were allowed. Event bursts were common in these observations, particularly those obtained 
in the March 2001 time period, so we carefully checked the cleaned lists to be sure that no obvious 
signal from these events remained after screening. 

Next, we combined the screened exposure lists for each channel within each observation using 
the default spectral registration provided by the software. This default assumption is necessary 
since there is not enough signal in the individual exposures to reliably cross-correlate the positions 
of narrow spectral features; the typical exposure time of ~ 2 — 3 ks per exposure yielded an average 
of ~ 2 — 3 counts per resolution element in the highest sensitivity (LiFl) channel and ~ 1 — 2 counts 
per resolution element in the SiC channels. The summed channel lists for each observation were 
corrected for geometric distortions, Doppler shifts, thermally-induced grating motions, astigmatism, 
detector backgrounds, and scattered light following the standard procedures (Sahnow et al. 2000). 
The extracted spectra for each channel in each observation were then flux calibrated and wavelength 
calibrated. We registered the spectra from each of the nine observations to a common wavelength 
scale by noting the integral pixel shifts required to align the same narrow absorption features 
observed in each spectrum. These shifts were typically less than 3 detector FUSE pixels and 
always less than 15 pixels. We then co-added the aligned spectra for all 9 observations to produce 
a final composite spectrum for each channel. 

We registered the composite channel spectra to a common heliocentric wavelength scale by 
performing a channel-to-channel registration similar to that performed at the observation level and 
then setting the zero point of the wavelength scale to that provided by the STIS spectra (discussed 
below). For the FUSE-to-STIS registration, we compared the velocities of lines of the same species 
observed by each instrument (e.g., C II AA1036.337, 1334.532; Si II AA1020.699, 1304.370; Fe II 
AA1144. 938, 1608.451). The C II and Fe II comparisons used lines with approximately the same 
line strengths, /A, to ensure that the comparisons were not influenced by asymmetric profiles that 
saturate in one line but not the other. Additional checks were performed using cross-element 
comparisons (e.g., Si II A1020.699 with S II AA1250.584, 1253.811, etc.). We then corrected the 
FUSE spectra to the Local Standard of Rest (LSR) reference frame. For the PG 1259-1-593 sight 
line, the correction for standard solar motion wlsr = ^hclio + 10.5 kms~^, assuming a solar speed 
of 19.5 kms-i in the direction l^ta = 56°, bstd = 23° (Delhayc 1965; Lang 1980; see also Kerr Sz 
Lynden-Bell 1986). This correction is within ~ 1.6 kms~^ of the LSR reduction based on a solar 
speed of 16.5 kms~^ in the direction / = 53°, 5 = 25° (Mihalas & Binney 1981). Unless stated 
otherwise, all velocities quoted in this paper are in the LSR reference frame. The FUSE spectra 
for PG 1259-1-593 have a nominal zero-point velocity uncertainty of ±5 kms~^ (Icr) after these 
calibrations have been applied. 

The FUSE data are oversampled in the spectral domain. Therefore, after calibration we binned 
the spectra to a spectral bin size of 4 pixels, or ~ 0.025 A (~ 7.5 kms"^). This binning provides 
approximately 3 samples per spectral resolution element of 22-25 kms~^. The night-only data 
used in this study have continuum signal-to-noise ratios S/N ~ 25 and 16 per spectral resolution 
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element at 1030 A in the LiFl and LiF2 channels, respectively, and S/N ~ 10 and 12 at 950 A in 
the SiCl and SiC2 channels, respectively. 

The version of CALFUSE (v2.2.2) used for this study has significant improvements over earlier 
versions used in previous studies, most notably improved wavelength solutions and background 
corrections tailored for day+night or night-only extractions. Portions of the fully reduced night- 
only SiC2 and SiCl spectra in the 915-955 A wavelength range are shown in Figure 1. The locations 
of the Lyman series H I and D I lines are indicated in the top panel along with the locations of 
prominent interstellar O I and N I lines. Terrestrial airglow emission features, which are present 
in all of the H I lines observed, are marked below the top spectrum with crossed circle symbols. 
Terrestrial O I and N I emissions are undetectable in the night-only data shown. The background 
subtraction for both channels results in nearly zero residual flux in the cores of the H I lines as 
expected, except at velocities affected by airglow contamination. (A similar conclusion holds for the 
very strong C II A 1036. 337 line observed in all four channels.) The slight rise in the sub-Lyman- 
limit regime (A < 912 A) in the SiCl channel may be the result of the radiative recombination 
continuum emission from atomic oxygen near 911 A (see Feldman et al. 1992; Lopez-Moreno et 
al. 2001) or improper background subtraction at these very short wavelengths. The SiC2 channel 
extends only to 916.6 A, so no independent check of the SiCl behavior at A < 916 A is possible. 

In the discussions that follow, we restrict our analysis to the orbital night-only data to minimize 
the impact of terrestrial H I and O I airglow emissions on the absorption lines of interest in this 
study. In particular, we concentrate on data from the two SiC channels covering the FUSE bandpass 
below 1000 A, but we also use the independent data from the two LiF channels to check the quality 
of the data at longer wavelengths. The checks allowed by multiple channel observations of the 
same spectral region are important for assessing noise in the data. When we derive H, D, and 
O abundances, we use both the SiCl and SiC2 data. Additional illustrations of FUSE spectra of 
PC 1259-F593 can be found in Richter et al. (2001b) and Wakker et al. (2003). 

2.2. Space Telescope Imaging Spectrograph Observations 

PG1259-I-593 was observed several times with HST/STIS as part of Guest Observer program 
GO-8695 on 2001 January 17-19 and 2001 December 19. An observation log is provided in Table 2. 
All observations used the intermediate-resolution echelle mode (E140M) and the 0^'2 x 0'.'06 slit 
to minimize the power in the wings of the spectral line-spread function (see Figure 13.90 in the 
STIS Instrument Handbook, Leitherer et al. 2002). This instrument mode provides a resolution 
of ~ 7 kms~^ (FWHM) per two-pixel resolution element and covers the 1150- 1729 A wavelength 
band with only five small gaps between echelle orders at A > 1634 A. For further details on the 
design and performance of STIS, see Woodgate et al. (1998) and Kimble et al. (1998). 

We reduced the data for each observation in the manner described by Tripp et al. (2001), 
including application of the two-dimensional scattered light correction developed by the STIS Team 
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(Landsman h Bowers 1997; Bowers et al. 1998). We weighted the extracted flux-caUbrated spectra 
for each observation by their inverse variances and added them together to produce a composite 
spectrum. The continuum signal-to-noise ratio per resolution element in the final co-added spectrum 
ranges from ~7 to 17 and peaks near 1400 A. We corrected the STIS data to the LSR reference 
frame by adding -1-10.5 kms~^ to the nominal heliocentric velocity scale provided by the standard 
processing. The zero-point accuracy of the STIS velocity scale for these data is approximately ±1 
kms~^. 

Examples of continuum normalized profiles for the absorption lines of several low ionization 
species are shown in Figure 2. Additional examples of absorption lines in these STIS spectra of 
PG 1259-F593 can be found in Richter et al. (2003) and Collins et al. (2003). 

2.3. H I 21 cm Observations 

High velocity clouds often show structure at angular scales down to at least ~l' (e.g., Wakker 
& Schwarz 1991), particularly in cloud cores. The PG 1259+593 sight line samples the Complex C 
core named CIII (Giovanelli, Verschuur, &; Cram 1973; Wakker 2001). The Leiden-Dwingeloo 
Survey (LDS, Hartmann k, Burton 1997) shows that this core has an extent of about 1° x 2° at a 
column density of 3 x 10^^ cm~^, with PG 1259-1-593 lying at its edge. To determine a more accurate 
estimate of the amount of high velocity H I in this region of the sky, we observed CIII with the 
Wcstcrbork Synthesis Radio Telescope (WSRT). Since this interferometer filters out the large-scale 
structure, we supplemented the WSRT data in the direction of PG 1259-1-593 with an Efi"elsberg 
single-dish observation having a 9f7 beam. This latter observation is described by Wakker et al. 
(2001). 

The half-power beam width (HPBW) of the WSRT primary beam is 35', which allows an area 
with a diameter of about 44' to be mapped. We therefore observed CIII with a mosaic of 2 x 4 
pointings spaced 27' apart, resulting in mapping over an area of (22'-|-27'-|-22') x (22'-|-3 x 27'-|-22') = 
71' X 125' arc minutes. To increase the sensitivity, we included an extra pointing in the direction of 

PG 1259-1-593 itself. For each of these pointings, we obtained full uv-planc coverage in 4 x 12 hours, 
using shortest spacings of 36, 54, 72, and 90 m. The correlator was set to cover the LSR velocity 
range between —333 and -M97 kms~^ with 2.1 kms~^ velocity resolution after on-line Hanning 
smoothing was applied. The WSRT observations were completed in late April 2001. 

R. Braun of the Netherlands Foundation for Research in Astronomy performed the calibration 
of the observations. The calibrated data were mapped using uniform weighting, and included a 
Gaussian taper in the uv-plane such that the final synthesized beam is 1' x 1'. Next, a map of 
the continuum was created for each pointing by adding the channels without H I emission. After 
subtracting the continuum, the individual pointings were cleaned using the Multi-Resolution 
Clean algorithm (Wakker & Schwarz 1988). To increase the accuracy of this step, the areas 
containing signal were first delineated, taking into account the overlap between pointings. The 
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resulting cleaned maps were then mosaiced together. The fully processed data have a final root- 
mean-square (RMS) residual of ~ 3.2 mJy beam~^, or ~ 0.5 K. 

Figure 3 shows a grey-scale map of N(H I) for Complex C core CIII based on the WSRT data 
(integrated between —148 and —109 kms^^). PG 1259+593 lies in the brightest concentration in 
the field, which coincides with the brightest spot in the LDS data. Although it seems unusual that 
the extragalactic background source lies toward the center of a cloud core, this is not an artifact 
of the mosaic process; if the pointing centered on PG 1259-1-593 is removed, the core is still clearly 
seen at the half-power point of each of the four surrounding pointings, and the final H I map looks 
similar, though noisier. We discuss the H I data and the derived column density for the high 
velocity gas in §6.1. 

3. The PG 1259-1-593 Sight Line 

PG 1259-F593 lies in the direction I = 121?09, b = +57?80 behind HVC Complex C, which 
spans Galactic longitudes between / ~ 30° and / ~ 150° in the northern Galactic hemisphere. Maps 
of the H I 21 cm emission in Complex C can be found in Wakker (2001) and Sembach et al. (2003); 
these maps show larger regions than the area around the sight line covered in Figure 3. Complex C 
has a mass M > 1.2 x 10^ Mq and a distance d > 5 kpc (or altitude z > 3.5 kpc) (van Woerden 
et al. 1999). Wakker (2001) suggests a distance limit d > 6.1 kpc and a mass M > 3 x 10^ Mq. 
Various determinations of the metallicity of the high velocity gas yield Z/Zq ~ 0.1 — 0.3 on a linear 
scale where Z = Zq is solar (Wakker et al. 1999; Gibson et al. 2001; Collins et al. 2003; Tripp ct 
al. 2003). Previously Richter et al. (2001b) found Z/Zq = 0.093±|]:J|f for Complex C based on a 
subset of the oxygen absorption lines toward PG 1259-1-593 considered here. 

In the direction of PG 1259-1-593, interstellar gas within 10 kpc of the Galactic plane that is 
participating in differential Galactic rotation has velocities vlsr ~ —30 to kms~^. There is also a 
large intermediate-velocity cloud known as the Intermediate- Velocity Arch in this general direction 
(Kuntz &; Danly 1996). Complex C has a velocity vlsr ~ —157 to —100 kms~^, which makes 
it possible to distinguish between the absorption produced by the HVC and the lower velocity 
foreground absorption produced by the nearby ISM {v^sR ~ ~5 kms~^) and the Intermediate 
Velocity Arch (IV Arch; vlsr ~ ~55 kms~^). The three principal groups of gas along the sight line 
are reasonably well separated in velocity in the H I 21 cm emission and weak ultraviolet absorption- 
line profiles shown in Figure 2. In the strongest absorption lines shown (Si II A1526.707, O I 
A1302.168, C II A1334.532), the ISM and IV Arch components blend together. It is important to 
note that the 21 cm emission associated with Complex C is stronger than that from the IV Arch or 
the Milky Way ISM in this direction, a point we will return to in the derivation of the Complex C 
H I Lyman series velocity structure in §6. Other sight lines through Complex C usually exhibit 
much weaker HVC 21cm emission and/or have a much more complicated velocity structure that 
hinders a clean separation of the Complex C gas from the ISM and intermediate velocity gas. For 
examples of the 21 cm emission toward other sight lines passing through Complex C, see Sembach 
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et al. (2003) and Wakker et al. (2001). The "enhanced" 21cm emission toward PG 1259+593 is 
undoubtedly due to the fact that the sight hne is located directly behind the Complex C core CIII 
(see §2.3). 

In addition to neutral gas, Complex C contains ionized gas that can be traced through its Ha 
emission (Tufte et al. 1998; Wakker et al. 1999), O VI absorption (Murphy ct al. 2000; Sembach 
et al. 2000, 2003), and C IV and Si IV absorption (Fox et al. 2003; see also §6 below). We must 
consider the impact of these ionized regions on our derivation of D/H in Complex C since they may 
contain trace amounts of H I detectable in the H I Lyman series lines (see §6) . 

4. General Methodology 
4.1. Overview 

Even though the PG 1259+593 sight line has many desirable attributes that facilitate a detec- 
tion of deuterium in Complex C, the conversion of this detection into a reliable estimate for the 
deuterium abundance depends upon many factors. It is necessary to take a methodical approach 
to determining the D/H and D/0 ratios in the high velocity gas. The sight line is more compli- 
cated than short sight lines through the local ISM used to determine the local value of D/H (e.g.. 
Moos ct al. 2002 and references therein; Hebrard & Moos 2003). It is also more complex than the 
high-rcdshift sight lines, which can be selected for simple intcrgalactic absorption velocity structure 
and typically contain one or two predominant absorption components (e.g., O'Meara et al. 2001). 
The velocity structure of the PG 1259+593 sight line does not allow for a precise measure of the 
H I column density in Complex C solely through analysis of the H I Lyman-series absorption due 
to blending of the high velocity H I absorption with the ISM and IV Arch components. However, 
it is possible to derive an accurate H I column density from 21 cm emission, provided that the 
angular resolution is sufficiently high to assure that small-scale structure within the radio beam 
is not biasing the result. It is also possible to draw upon the information in the O I absorption 
profiles for the sight line and to constrain the range of possible Lyman-series absorption parameters 
for both D I and H I. 

Our adopted method for determining the D/H and D/0 ratios in Complex C can be broken 
down into several key steps: 1) determination of an accurate O I column density for Complex C 
and a model for the velocity structure of the ncTitral gas along the PG 1259+593 sight line using 
the numerous O I absorption lines in the FUSE and STIS spectra, 2) application of the O I velocity 
model to the H I Lyman-series lines in the FUSE data to create an H I model that reproduces the 
observed H I absorption and 21 cm emission profiles, 3) determination of the H I column density 
in Complex C from the interferometric H I 21 cm data described in §2.3, 4) determination of the 
amount of D I required to account for the Complex C absorption in the negative velocity wings of 
the H I Lyman-series lines, and 5) consideration of the uncertainties associated with steps (1) - (4). 
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4.2. Methods 

We use a variety of spectral line analysis techniques to determine the velocity structure and 
column density of each species studied. Our primary tool for quantifying the properties of the 
absorption lines in the FUSE and HST/STIS spectra is a suite of software for fitting line profiles, 
written specifically for this study in the Interactive Data Language (IDL). The software constructs 
synthetic line profiles that can be compared to the observed data after convolution with an instru- 
mental line spread function. It allows complex absorption lines to be modeled as a superposition 
of components with Maxwellian velocity distributions, each of which can be described as Voigt 
functions with appropriate natural damping constants. Each component has a central velocity 
{{vi)), line width (bj), and column density (Nj), with the best-fit parameters found through a x^- 
minimization of the differences between the model profiles and the data. To explore the sensitivity 
of our adopted results to various input parameters/assumptions and boundary conditions, we per- 
form some of these profile fitting analyses multiple times with other profile fitting codes used in 
FUSE analyses of other sight lines (e.g., the Owens. f code - see Hebrard et al 2002; Lemoine et al. 
2002). In all cases checked, the column densities are consistent to within the la errors, confirming 
that the formal uncertainties have been estimated properly. For the O I and D I lines we also use 
curvc-of-growth analyses to estimate column densities and to check the reliability of the profile 
fitting results. We describe these different analyses in the following discussions of the absorption 
produced by each species (H I, D I, O I) . 



4.3. Input Data and Reference Abundances 

The input atomic data parameters for this study are known with high enough accuracy that 
they do not present a significant source of systematic uncertainty in our final results. We adopt 
wavelengths, oscillator strengths (/-values), and radiation damping constants from the atomic data 
compilations of Morton (1991, 2003). For O I, the original sources of the /-values are Zeippen, 
Seaton, & Morton (1977), Biemont & Zeippen (1992), and Tachiev & Froese Fischer (2002). The 
source of /-values for H I and D I is Pal'chikov (1998). Wc use multiple transitions of H I, D I, and 
O I in our line analyses. Thus, typical uncertainties of 10 — 20% in the /-values for individual lines 
are mitigated by considering several transitions simultaneously in the column density determina- 
tions. To identify molecular hydrogen lines and assess their possible contamination of the atomic 
absorption features, we use the H2 line lists of Abgrall et al. (1993a,b). 

Throughout this work, we adopt a solar reference abundance (O/H)© = 4.90 x 10~^ from 
AUende Prieto, Lambert, Sz Asplund (2001). This reference solar oxygen abundance is in good 
agreement with the value of 4.5 x 10~^ derived from analyses of the solar [O I], O I, and OH 
line shapes and asymmetries implied by new 3-D hydrodynamical model atmospheres (Asplund 
2003). It is in better agreement with the nearby average ISM oxygen abundance of (0/H)ism = 
3.4 X 10~^ (Meyer 2001; Meyer, Jura, & Cardelli 1998) than earlier determinations based on 1-D 
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model atmospheres. It also agrees well with recent determinations of H II region oxygen abundances 
(see Pilyugin, Fcrrini, & Shkvarun 2003 and references therein). We will adopt this solar reference 
abundance in our discussions of the metallicitics of the various absorbers along the sight line. Using 
(0/H)q = 5.45x10"^ from Holweger (2001) or (O/H)© = 7.41x10"^ from Grevesse & Noels (1993) 
yields values of [0/H] that are a factor of 0.05 dex and 0.18 dex lower, respectively, than those 
adopted in this study. 

4.4. Simplifying Assumptions 

4.4- 1- Number of Components in the Model 

Throughout this work, we adopt the minimum number of velocity components required to 
produce acceptable fits to the observed absorption lines in the FUSE and HST/STIS observations. 
Our choice of component structure is guided in part by the high-resolution HST/STIS data available 
for several metal-line species (e.g., O I, Si II, S II, Fe II). In all, 8 components are required to fit 
the O I, H I, and D I profiles. The attributes of these components will be discussed for each species 
in §§5-7. A standard F-test indicates that adding additional components does not significantly 
improve the quality of the fit to the observed data, with the possible exception of an additional 
weak H I absorption feature discussed in §7.1. Adopting fewer components does not provide an 
acceptable fit to both the O I and H I absorption-line data. Some of the identified components, 
especially those at low velocities, are likely to consist of multiple sub-components that cannot be 
resolved at the resolution of the FUSE and HST/STIS data. 

4-4-2- Velocity Structure Assumptions in the Model 

We require the O I, H I, and D I lines to have a similar velocity structure, but we do not 
simultaneously fit these species with other metal-line species (e.g., N I, Si II, S II, Fc II) as has 
been done in some previous D/H analyses of other sight lines. To do so would impose additional 
constraints on the O I, H I, and D I lines that are not only unwarranted but also potentially 
misleading for this sight line. The absorption profiles of other species may be influenced significantly 
by the presence of ionized gas in the different components along the sight line. Unlike sight lines 
to nearby stars where the coupling of many species makes sense because the velocity structure of 
the ISM is simple and the neutral and ionized species have very similar velocity structure (e.g., 
Vidal-Madjar et al. 1998; Hebrard et al. 2002), the PG 1259-^593 sight fine contains a variety of 
absorbing regions with a range of physical properties. Linking the O I velocity structure directly 
to the profile shapes of species that trace both neutral and ionized gas could result in an inferred 
velocity extent of the neutral Complex C absorption that is sensitive to the profile shapes of singly 
ionized species in Complex C or in the IV Arch gas at nearby velocities. A secondary concern with 
tying the velocity structure of O I to ions such as Si II or Fe II is the inclusion of refractory elements 
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into dust. Differential depletion of elements into dust grains or selective shock-disruption of grains 
can alter the relative gas-phase abundances of ISM clouds encountered along the sight line, which 
in turn can translate into different absorption-line velocity structures for different species. We use 
all of the information available to construct an initial guess at the velocity structure of the sight 
line, but take the conservative approach of coupling the velocity structure of only those species 
(i.e., H I, D I, O I) for which ionization and dust depletion effects are unlikely to differ significantly. 
Information for the singly ionized species can be found in recent studies by Richter et al. (2001b) 
and Collins et al. (2003). 

4.4.3. The FUSE Line Spread Function 

Wc assume that the FUSE line spread function (LSF) has a Gaussian shape and a constant 
width (FWHM) of 0.075 A 12 FUSE detector pixels) throughout the FUSE bandpass. This 
width corresponds to a velocity width of FWHM 20 — 25 kms~^, depending upon wavelength. 
This is slightly broader than the nominal value of 20 kms~^ often assumed in FUSE studies (see 
Sahnow et al. 2000). It is reasonable to expect the LSF for the combined PG 1259-1-593 data 
to be degraded slightly compared to the nominal LSF for single FUSE exposures because the 
number of individual exposures in the combined set of observations is large and registration of the 
individual exposures is certain to introduce some minor spectral degradation, particularly in the 
SiC channels. Given the velocity structure of the sight line, it is difficult to find narrow isolated lines 
to determine the LSF width directly. The narrowest lines for which this is possible in the FUSE 
spectrum of PG 1259-1-593 have observed (instrumentally convolved) widths of 0.08 A (LiFlA 
detector segment) and ~ 0.10 A (SiC2A detector segment). The measured widths are broader than 
the adopted LSF width but are consistent with the expected line widths based on the velocity 
structure observable at higher resolution in the STIS data shown in Figure 2. 

In previous investigations of D/H along kinematically simple sight lines, a multiple-component 
LSF has sometimes been used to fit the data. For example, Kruk et al. (2002) and Wood et 
al. (2002) used a two-component LSF with a narrow component (FWHM = 9 pixels) containing 
~ 70 — 75% of the total LSF area and a broad component (FWHM = 17 — 24 pixels) containing 
~ 25 — 30% of the LSF area. When summed, the double-Gaussian functions studied have a roughly 
Gaussian shape with a width similar to that adopted in this study. The resulting column densities 
of D I and O I obtained with a double-component LSF have been similar to those obtained with 
single-component LSFs (Kruk et al. 2002; Wood et al. 2002). The reasons for this are straight- 
forward: the nearby absorption is simple, and the effects of a broad LSF component are most 
pronounced for strongly saturated lines requiring very accurate estimates of the optical depth at 
low residual flux levels. Toward PG 1259-1-593, the velocity structure is much more complex. Since 
we must consider the case where weak absorption is in close proximity to strong absorption, it is 
not so easy to dismiss the effects of a non-Gaussian LSF. For our PG 1259+593 FUSE datasct, 
which is the combination of many individual exposures, the central-limit theorem implies that the 
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LSF of the combined data should also closely represent a single Gaussian function. After modeling 
the sight line and considering the impact of a more complex LSF, we find that the addition of a 
second LSF component with broad wings docs not change our estimates of the H I, D I, and O I 
line strengths. Other systematic uncertainties dominate the potential uncertainties caused by the 
adoption of a single-component LSF rather than a two-component LSF of the type used in previous 
FUSE studies. 

We performed a consistency check on our LSF estimation by comparing the shape of the 
Fe II A1144.938 line in the FUSE LiF2A spectrum with that of the Fe II A1608.451 line in the 
HST/STIS spectrum. The two lines have intrinsic strengths (/A) within ~ 20% of each other, so 
direct comparisons of the velocity structure and component absorption depths could be made by 
relating the FUSE spectrum to various smoothed versions of the higher resolution STIS spectrum. 
We found that a good approximation to the FUSE Fe II profile was obtained by convolving the 
STIS spectrum with a Gaussian function having FWHM = 20±4 kms~^. (At 1145 A, our adopted 
instrumental width of 0.075 A corresponds to FWHM = 19.6 kms~^.) Smaller smoothing widths 
left structure in the smoothed STIS spectrum that was not obvious in the FUSE data, while larger 
smoothing widths underestimated the line depths and overestimated the line widths of the FUSE 
absorption. This comparison docs not reveal the exact shape of the FUSE LSF, but it is sufficient to 
bracket the LSF width and demonstrate that the adopted Gaussian LSF width is well-constrained. 



4 -4 ■4- Estimation of the Quasar Continuum 

Another simplifying assumption we make is that there are no undulations in the ultraviolet 
continuum of PG 1259-1-593 on scales of less than a few Angstroms. This assumption is made in 
essentially all studies of quasar absorption line spectra. Compared to the spectra of normal and 
degenerate hot stars, power-law quasar continua are exceedingly smooth over small wavelength 
intervals; see, for example, the FUSE spectrum of the quasar 3C 273 (Sembach et al. 2001b). The 
PG 1259+593 spectrum contains no signatures of a stellar population contribution to the continuum 
light, which can be difficult to model for some AGNs, particularly Type II Seyferts. We model the 
PG 1259-1-593 continuum in the vicinity of the interstellar lines as a smoothly varying function of 
wavelength. The low-order Legendre polynomial fits adopted for the SiC channel data below 960 A 
are shown in Figure 1. The quasar continuum is relatively flat (i.e., dFx/dX ~ 0) over the large 
wavelength interval considered (~ 900 — 1100 A). Uncertainties associated with the placement of 
these continua are small compared to the other sources of uncertainty discussed below. Continuum 
placement is especially important for the weak lines of D I, for which we have estimated upper 
limits including this systematic uncertainty. We used the methods outlined by Sembach & Savage 
(1992) to determine the continuum levels and to extract the continuum placement uncertainties. 
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4.4-5. Possible Contamination from Absorption Lines of Other Species 

There is little molecular hydrogen along the PG 1259+593 sight line, as evidenced by the 
paucity of H2 R(0) and R(l) Lyman-band absorption features in the FUSE spectra. Richter et al. 
(2001a) estimate Ej=o logNj(H2) < 13.96 {3a) in Complex C. There may be a small amount of 
H2 absorption present in the strongest lines of the J = — 3 rotational levels at the velocities of 
the Milky Way ISM (i^lsr ~ —5 kms^""^); we estimate E^=o ^°S^j(-'^2) ~ 14.7. A small amount 
of H2 may also be present in the IV Arch [ Ej=o ^j(H2) = 14.0ibQ^ij ]. We examined each O I 
and D I line used in this study for possible interstellar H2 lines that might be confused with the 
atomic absorption in Complex C at —160 < vlsr < — 110 kms~^. We found that the predicted 
interstellar H2 lines in the J = — 3 rotational levels with the total column densities quoted above 
produce negligible absorption in the vicinity of these lines. H2 absorption by the foreground ISM 
can safely be ignored in our analysis of the Complex C absorption features. Some II2 lines arising 
in the Milky Way ISM occur at wavelengths corresponding to the low-velocity portions of the H I 
lines, but these features do not affect the inferred strengths or shapes of the strongly saturated H I 
lines. 

There are no known absorption lines produced by the intergalactic medium (IGM) along the 
PG 1259+593 sight line at the wavelengths of the O I and D I lines between 915 A and 950 A. The 
strongest IGM systems occur at redshifts z = 0.04606 and 0.21949 (Richter et al. 2003). In both 
cases the H I Lyman-series lines in these systems occur at wavelengths longward of the 910 — 950 A 
spectral region. The remaining 19 IGM systems with z = 0.00229 — 0.43148 have H I detectable 
only in the Lya — Ly7 transitions (Richter et al. 2003), so there is no possibility of contamination 
at shorter wavelengths by those systems. Of the possible metal-line species that are likely to be 
detectable, C III A977.020 tends to be the strongest and most common; it too is safely longward 
of the absorption features considered in this study. Potential redshifted extreme ultraviolet species 
(e.g., O II A834.466, O III A832.927, O IV A787.711) in the intergalactic absorbers along the sight 
line do not pose a problem either. 

5. Oxygen 

5.1. 01 Velocity Structure 

The first step in our study of the D/H ratio in Complex C is a determination of the O I velocity 
structure of the gas along the entire sight line in the velocity range luisijl < 300 kms~^. We use 
both the moderate-resolution FUSE observations of numerous O I lines at far-ultraviolet wave- 
lengths, combined with the higher resolution HST/STIS echelle observations of the O I A 1302. 168 
line, to examine the optical depth of the absorption in detail. (Note that the weak intersystem O I 
line at 1355.598 A is not detected in the STIS spectrum.) Deriving the O I velocity structure is 
critical because it provides an excellent initial approximation to the velocity structure of the D I 
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and H I lines and is required for an accurate metallicity ([0/H]) estimate for the Complex C gas. 

The ionization fraction of O I is coupled to that of H I (and thus D I) through a strong charge 
exchange reaction at typical interstellar densities (Field & Steigman 1971): 

0+(^S3/2) + H° ^ 0°(=^P2) + H+ + 0.02 eV. (1) 

Detailed CLOUDY photoionization models for a range of physical conditions spanning those likely to 
be encountered along the PG 1259+593 sight line show that the ionization fractions of O I and H I 
track each other closely. For the specific case of the 3C351 sight line through Complex C, Tripp et 
al. (2003) find that the ratio of ionization fractions of H I and O I is unity to within a few percent 
for —5.0 < \ogU < —3.0, where U = n-y/nn is the ratio of hydrogen-ionizing photon density to 
total hydrogen number density. A valuable check that the ionization balance of O I and H I is not 
impacted by extreme conditions is provided by the observed ratio of Ar I to O I in Complex C 
toward PG 1259+593. Collins et al. (2003) find that the two species have a ratio of column densities 
within a factor of a few of the expected value of Ar/0 in a solar-abundance gas. Neither Ar nor O 
is severely depleted into dust grains, and since both are a-elements, their chemical histories should 
be similar. Tripp ct al. (2003) deduce that logC/ ~ —4.0 based on the Collins et al. result (see 
their Figure 14). Therefore, the ionization correction for O I should be small. 

The present-day cosmic abundances of D and O are within a factor of ~ 20 — 50 of each other, 
which makes the weak O I lines in the FUSE bandpass an excellent guide to the velocity structure 
expected for the D I lines. The ratio of N/A for the O I lines to the D I lines considered in this 
study is roughly unity. The strong O I A1302.168 line in the HST/STIS data helps to constrain 
the amount and velocity structure of neutral gas contained in low column density components 
unobservable in the weaker O I transitions in the FUSE bandpass, which is an especially important 
consideration in modeling the H I velocity structure at the velocities of the D I lines. An additional 
incentive for determining the O I velocity structure is that some of the O I lines overlap the D I and 
H I lines (see Figure 1). It is therefore desirable to estimate the contributions of the O I absorption 
to the blended (D I+O I) absorption features observed toward PG 1259+593. 

We list the O I lines considered in this study in Table 3. The table entries include the 
wavelength, line strength (log /A), instrument used for the measurements, and notes about the 
line and possible blends with other absorption features. We have adopted /-values from Morton 
(1991, 2003), but using the /-values recommended by Verner, Barthel, &; Tytler (1994) does not 
change the predicted O I absorption strengths appreciably. The transitions considered span a factor 
of 276 in line strength (/A) from A922.200 to A1302.168. The lines used in our O I analyses below 
are indicated in column 5 of Table 3. 

To reproduce the observed velocity structure of the O I lines, we start with a model that 
contains components clearly present in the neutral gas. The HST/STIS data shown in Figure 2 
are the best guide available. Three primary components at vlsr = -129 kms~^ (Complex C), -55 
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kms-^ (IV Arch), and -2 kms~^ (Milky Way ISM) are prominent in the S II AA1250.584, 1253.811, 
1259.519 Hnes; these primary components are also seen in the H I 21 cm emission profiles shown 
at the top of the figure. A second Complex C component at -112 kms^^ is present in the O I 
A1302.168 line profile and in the strong Si II, Al II, and Fe II lines observed by STIS. An additional 
ISM component at +21 kms~^ is also visible in these same lines, and a weak intermediate velocity 
component is present at +69 kms^^ in several lines (C II A1334.532, O I A1302.168, Al II A1670.787, 
Si II A1260.422, Si III A1206.500). These 6 components form the basis for our initial model - a set 
of velocities, line widths, and estimated optical depths. The resulting model was convolved with 
the FUSE and STIS LSFs, and the parameters were varied to minimize the residuals between the 
model fits and the observed FUSE and STIS O I data. During this process, we found that two 
additional components at -82 and -29 kms~^ were needed to produce an acceptable fit to the 
ensemble of O I lines. These components are included in the final model spectrum, which has the 
parameters listed in Table 4. 

We show the results of the O I fitting process in Figure 4 for a representative set of O I lines 
with strengths -1.0 < log /A < 1.826. The STIS data for the A1302.168 line are shown at the 
bottom of the figure, and the FUSE data (orbital night-only) arc shown as either solid (SiC2) or 
dashed (SiCl) lines. The adopted model, which is overplotted on the data as a heavy solid line, 
reproduces the observed absorption features in both the weak and strong lines quite well. The 
predicted strengths of those lines that are not shown are equally well represented. This model 
solution is not necessarily unique, but it provides some fundamental information. First, it indicates 
that the Complex C absorption has two distinct components. Earlier determinations of N(0 I) using 
subsets of the present data (Richter et al. 2001b; Collins et al. 2003) assumed single-component 
representations for the absorption. Second, the Complex C absorption is distinct from the IV Arch 
absorption. Therefore, degeneracies in the parameters for the lower velocity gas are not critical for 
estimating N(0 I), although they are important for understanding the possible range of H I (and 
thus D I) velocity structure. Third, the model reproduces the absorption in both weak and strong 
O I transitions, so it provides a reliable estimate of the O I absorption for those remaining O I lines 
that are blended with other absorption features (see Table 3). 

In constructing the model shown in Figure 4, wc found that the reduced absorption strength 
a-t vlsr = 20.5 kms"^ {vheUo = 10 kms~^) in the HST/STIS O I A1302.168 profile must be due 
to geocoronal O I airglow emission. This feature is indicated by the cross-hatched region in the 
figure. The velocity of the emission feature matches (to within 1 kms~^) that of the geocoronal 
H I Lja in the same STIS spectrum. Additional weak airglow features of O I* A1304.858 and O I** 
A1306.029 at vlsr = 20.5 kms~^ are also present in the STIS spectrum at the l-2cr confidence level, 
lending further support to the geocoronal identification of the feature. The FUSE O I lines, which 
are weaker than A1302.168 and were observed during orbital night, are not expected to contain 
geocoronal emission. They show considerable absorption at the velocity of the emission feature, 
as do other metal lines in the STIS spectrum that are not substantial constituents of the Earth's 
atmosphere (e.g.. Si II AA1304.370, 1526.707; Fe II A1608.451). Thus, our O I model requires an 
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interstellar absorption component at this velocity. The precise details of how much interstellar O I 
absorption occurs at the velocity of the atmospheric O I A1302.168 airglow feature do not affect 
the conclusions of this paper. 

P II A1301.874 (log /A = 1.351; Morton 1991) occurs at -68 kms"^ with respect to O I 
A1302.168. Estimates of the line strengths of interstellar ~ kms"^) P II A1152.818 (log /A = 
2.435, Wx = 97± 16 mA) and P II A1532.533 (log /A = 1.067, Wx < 50 mA [3a]) imply a column 
density log N(P II) ^ 13.6 and a corresponding equivalent width of only « 10 mA for the 1301.874 A 
line. Therefore, the interstellar P II line is buried in the strong O I absorption shown in Figure 4 
and contributes negligibly to the observed O I profile. The same is true of the IV Arch absorption. 
P II A1301.874 absorption associated with the IV Arch would occur near an O I A1302.168 velocity 
of -123 kms^^. This line is also expected to be weak (< 10 mA) since the strength of the IV Arch 
P II A1152.818 absorption is comparable to that of the ISM line. As a consistency check, the 
IV Arch P II A1301.874 absorption should have a strength similar to that of P II A1301.874 in 
Complex C (Richter et al. 2001b), which has a velocity of ^ -197 kms"^ on the O I A1302.168 
velocity scale shown in Figure 4. There is no detectable absorption at this velocity over a 50 kms~^ 
wide interval (Wx < 21 mA, log N < 13.91 [3a]). We conclude that the P II A1301.874 lines for the 
ISM, the IV Arch, and Complex C do not contribute significantly to the observed O I A1302.168 
absorption shown in Figure 4. 

5.2. O I Column Densities 

In Table 4 we list the velocities, line widths, and approximate O I column densities for each of 
the components used in the fit to the O I lines. The errors on the column densities and widths are 
difficult to estimate precisely for some of the individual ISM components because of the complex 
velocity structure of the overlapping features. However, the summed column densities for the 
three primary absorption groups (Complex C, IV Arch, ISM) and the positive velocity IVC are 
well-constrained by the observations. The total column density and la error for each group are 
listed after the individual component values in Table 4. We discuss the column densities in each 
component group below. 

5.2.1. Complex C 

Most of the O I in Complex C is contained in the component at -129.5 kms~^ (component 
#1). This component dominates the absorption produced by Complex C in the weak O I lines 
in the FUSE bandpass and is the main constituent of the A1302.168 absorption observed with 
HST/STIS. The weaker component at -112.5 kms^^ (component ^2) accounts for the absorption 
in the lower velocity Complex C absorption in the A1302.168 line and A1039.230 lines, but does not 
contribute significantly to the weaker O I absorption features observed. The primary constraints 
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on the weaker feature are set by the shape of the A1302.168 hne since the absorption produces only 
a small inflection in the A1039.230 line. The widths of both Complex C components are tightly 
constrained by the steep absorption walls of the A1302.168 line. 

The large range of O I line strengths available constrains the total column density of the 
Complex C absorption. We find a total Complex C O I column density of 



where the errors are dominated by systematic effects associated with the saturation correction 
necessary for the higher velocity component containing most of the column density. The statistical 
error is negligible (< 5%) compared to the systematic error since many lines are used to determine 
the quality of the fit. The level of inferred saturation in the strong component must be judged 
primarily by the strengths of the lines observed in the FUSE bandpass because the A1302.168 
line has a very high optical depth. Fortunately, the width of the stronger component is well- 
constrained (b = 6.0 ±1.0 kms~^). The small amount of O I in the lower velocity Complex C 
component contributes at most ~ 10% to the total O I column density in Complex C. However, this 
component must be included in the profile fits because it affects the inferred b-valuc for the strong 
component; single-component fits to the ensemble of O I lines shown in Figure 4 systematically 
underestimate the true column density by ~ 20%. 

Previous studies of the O I abundance in Complex C used single-component curves of growth 
to estimate the O I column density and found N(0 I) ~ 5.7 x 10^^ cm~^, which is within our 
1(7 confidence range but about 20% less than our preferred value (Richter et al. 2001b: log N = 
15.77±g:|J; Collins et al. 2003: log N = 15.75±g:^|). Our best-fit column density estimate, log 
N(0 I) = 15.86 ± 0.15, is slightly higher than the values reported in these previous studies because 
we derive a smaller b-value for the main Complex C component than the single-component COG 
studies; statistically, all three values are within the quoted errors on each measurement. The two 
previous studies adopted b « 10 kms~^ for the Complex C absorption, compared to b 6 kms~^ 
for the stronger component in this study, and neither made use of the detailed velocity structure 
revealed by the high-resolution STIS data even though the A1302.168 line was included in the COG 
fits. One expects a larger b-value for the single-component COG fits because the weaker (—112.5 
kms~^) Complex C component contributes to the total equivalent widths of the strong lines used 
in the COG fits (AA1302. 168, 1039.230). In these cases, the b-value is an "effective" b-valuc for 
the combination of the two components weighted by the individual line strengths, line widths, and 
velocity separation of the two components. The Complex C absorption requires a minimum of two 
components to fit the non-symmetrical velocity structure observed in the A1302.168 line. Adding 
additional weak components does not improve the fit in a statistically meaningful way. 

We checked the veracity of our O I profile fit results by measuring the equivalent widths of 
the lines shown in Figure 4 and plotting these data points on curves of growth. The measurements 



N(0I) = (7.2 ±2.1) X 10^5 
N(0I) = (7.2 ±4.2) X 10^^ 



cm ^ (68% confidence) 
cm~^ (95% confidence) 
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are listed in Table 5, and the curves of growth are shown in Figure 5. We considered two cases: a 
single-component curve of growth (like the ones used in previous analyses) and a double-component 
curve of growth based on the fit parameters listed in Table 4. For both curves of growth, wc have 
plotted the corresponding O I A1302.168 model absorption profile over the observed data shown in 
the inset box in each panel of Figure 5. The single-component fit yields log N(0 I) = 15.66±Q Qg 
and b = 9.7 it 0.7 kms""^. This column density is at the 2(7 lower limit of the column density in 
our best fit model. The double-component curve of growth improves upon the single-component 
column density estimate by providing both a good fit to the observed equivalent widths and a 
better fit to the velocity structure of the A1302.168 line profile, as can be seen from the insets in 
Figure 5. 

The error derived for N(0 I) from the profile fitting process is shown graphically on the double- 
component curve of growth illustrated in Figure 6. The dashed curves represent the COG results 
when the value of N(0 I) in the stronger Complex C component is changed by 4.2 x 10^^ cm~^ 
(±2(7). The change in column density affects the predicted strengths of the weak O I lines to a 
much greater extent than the strength of the A1302.168 line. The range of column density values 
allowed is consistent with the range predicted from the profile fits. 



5.2.2. IV Arch and Galactic ISM 

The column densities derived for the IV Arch and the ISM components toward PG 1259-1-593 
do not strongly affect the value of N(0 I) derived for Complex C. The total O I column density 
for the IV Arch in our model is log N(0 I) = 16.07 it 0.10, which is roughly a factor of 2 smaller 
than the value of 16.34ibQ;27 derived by Richter et al. (2001b) from a single-component COG. The 
two estimates are within the measurement uncertainties of each other, and the exact value adopted 
depends weakly upon the lower velocity integration limit. The two weak IV Arch components at 
-82 and -29 kms~^ are not visible in the H I 21cm emission profiles; the column densities are 
an order of magnitude or more lower than the column density in the principal component at -55 
kms~^. However, these two components contribute to the O I (and presumably H I and D 1) 
absorption. Using log N(H 1) = 19.48 derived from the Effelsberg data shown at the top of Figure 2 
(see also Wakker et al. 2003), we find [0/H]ivArch = log [N(0 I)/N(H I)] - log (O/H)0 ^ -0.10. 

The total O I column density derived for the ISM components with vlsr = -2.5 and -|-21 
kms~^ is log N = 16.11 ± 0.10, which is in reasonable agreement with what would be expected 
for gas with solar abundances and an H I column density log N(II I) = 19.67 derived from the 
Effelsberg data shown in Figure 2. We find (0/H)ism ~ 2.75 x 10"^, or [0/H]ism ~ -0.25, for the 
gas between -30 and -|-30 kms~^. 
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6. Hydrogen 

6.1. H I 21 cm Emission Data and the H I Column Density of Complex C 

To derive the best value for N(H I) in the direetion of PG 1259+593, we eombined the WSRT 
map shown in Figure 3 with the single-dish Effelsberg observation centered on the QSO (see §2.3). 
The method we followed is fully described by Schwarz & Wakker (1991) and Wakker, Oosterloo, &: 
Putman (2002). In summary, we first "observed" the WSRT map with a 9'.7 beam. The resulting 
spectrum was subtracted from the Effelsberg spectrum, and the full-resolution WSRT data were 
added back in. To each of these spectra we fit a single Gaussian emission feature to estimate the 
H I column density in Complex C. Figure 7 shows the results of these procedures. From the fit 
shown in the figure, we estimate N(H I) = (8.9 it 0.7) x 10"*^^ cm~^ in Complex C in the direction of 
PG 1259+593. The single-component fit to the Complex C emission has FWHM = 16 + 2 kms~^ 
(b = 9.6 + 1.2 kms~^). A straight integration of the H I profile from -145 to -109 kms"^ yields 
N(H I) = (9.0 + 0.6) X 10^^ cm"^, where the statistical error in this case is given by propagating the 
RMS noise per channel (0.32 Jy beam~^) over 28 channels of width 1.288 kms~^ per channel. Note 
that N(H I) is nearly identical to that derived from the Effelsberg data because N(H l)^^,,^^^^^^ ~ 
N(H I) interferometer- ^hc higher rcsolutiou WSRT data reveal that the profile is about 20% narrower 
than previously measured, showing that some of the width of the single-dish profile is due to a 
small velocity gradient in the 9'.7 Effelsberg single-dish beam. 

Sources of systematic errors associated with our determination of N(H I) in Complex C include 
cleaning and calibration of the data, scaling of the single-dish data, the shape of the single-dish 
beam, the "negative bowl" in the interferometer maps, and the fine-scale H I structure within the 
beam. Adding the uncertainties associated with the first four of these errors together yields an 
"instrumental" systematic error of 0.7 x 10^^ cm~^ if all act in the same sense; this uncertainty is 
comparable to the statistical error quoted above. Although the WSRT data are not able to provide 
information about gas structures on scales smaller than the 1' effective beam size (~ 1.5 pc at a 
distance of 5 kpc), it is possible to assess the magnitude of the uncertainties caused by H I fine 
structure. We calculated the fluctuations in the H I column density in box sizes of 3 x 3 pixels 
(60' x 60'), 5x5 pixels (100' x 110'), and 7x7 pixels (140' x 154'). We find averages and variations 
of (9.0 + 0.5) X lO^'^ cm-2, (10.0 + 1.6) x 10^^ cm-2, and (11.0 + 2.4) x lO^^ cm^^ respectively. 
The minimum and maximum values in the 3x3 pixel box are 8.4 x 10^^ cm~^ and 10.0 x 10^^ 
cm~^, respectively. Some of these small-scale structure effects are already covered indirectly in 
the "instrumental" systematic error estimate, but we note that there could still be additional 
uncertainties contributed by fluctuations at small scales that we cannot assess. A total systematic 
error of ^ 1.1 x 10^^ cm~^ accounts for most of the systematic errors listed above, including small- 
scale structure within the 3x3 pixel box, even if the majority of these uncertainties affect the column 
density measurement in the same sense. We adopt this value as our 95% confidence uncertainty on 
the systematic error, and a value of 0.8 x 10^^ cm^^ as our 68% confidence uncertainty. 

Adding the corresponding systematic and statistical errors in quadrature to produce a total 
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error yields a final H I column density estimate based on the 21 cm emission data of 

N(H I) = (9.0 ± 1.0) X 10^9 cm-2 (68% confidence) 
N(H I) = (9.0 ± 1.6) X 10^9 cm-2 (95% confidence). 

We adopt this column density and its associated errors in our discussion of the H I Lyman-series 
profiles, and we use the 21 cm data as a constraint on the velocity structure of the neutral gas along 
the sight line. 

Other sight lines show higher velocity 21 cm emission in Complex C than the PG 1259+593 
sight line (e.g., the Mrk 876 sight line; Murphy et al. 2000). Therefore, it is useful to consider 
whether the PG 1259+593 sight line might also contain some lower column density H I at higher 
negative velocities that might be confused with D I in the Lyman-series lines. The 21 cm emission 
data set a fimit of N(H I) < 2.0 x 10^^ cm^^ at velocities —275 < vlsr < —175 kms^^, the velocity 
range over which low column density H I could masquerade as D I absorption. In constructing this 
limit, we have used the high-S/N side-lobe-cleaned 21 cm spectra shown in Figure 2 from the NRAO 
140-foot telescope (Murphy, Sembach, & Lockman, private communication). These data set a more 
stringent constraint on the amount of very high velocity gas than either the WSRT or Effelsberg 
data. 

A limit on the amount of H I that could be present at very high velocities can also be set by 
the lack of absorption in the O I A1302.168 and C II A1334.532 profiles. We find 3a equivalent 
width limits of 30 mA for both lines over this 100 kms~^ velocity interval. These limits correspond 
to N(0 I) < 4.1 x 10^^ cm~^ and N(C II) < 1.5 x 10^^ cm~^, assuming the lines lie on the linear 
part of the curve of growth. For gas with solar abundances, (O/H)© = 4.90 x 10""^ (Allende Prieto 
et al. 2001) and (C/H)e = 2.45 x 10"^ (Allende Prieto, Lambert, & Asplund 2002), the metal-line 
limits translate into H I limits of < 8.4 x 10^^{Z/Zq)-'^ cm'^ (from O I) and < 6.1 x 10^^{Z/Zq)-'^ 
cm^^ (from C II). For Z/Zq '-^0.1, these limits are only a factor of a few times more stringent 
than those provided by the H I 21 cm emission data. Consideration of narrower velocity intervals 
would produce smaller upper limits. Unfortunately, neither the 21 cm emission nor the metal-line 
absorption provide strict enough limits to rule out the possibility that some of the absorption 
attributed to D I may be high- velocity, low-density H I (see §7.1). 

While it is not possible to distinguish between D I and interloping H I in the Lyman-series 
lines observed by FUSE, an estimate of the amount of neutral gas at high positive velocities can 
be made. The lack of H I absorption at +100 < ^lsr < +300 kms^^ in the Lye line sets a limit 
of log N(H I) < 14.55 (3cr) on high positive velocity gas along this sight line. This is the strongest 
Lyman-series line for which this velocity interval is not confused by overlapping absorption features; 
Lya and Ly/3 have damping wings from lower velocity components, Ly7 has a nearby intervening 
intergalactic absorption line at these velocities, and LyS has nearby interstellar O I absorption at 
these velocities. 
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6.2. Lyman-Series H I Velocity Structure 

An important step in our estimation of the deuterium abundance and the D /H ratio in Com- 
plex C is the conversion of the O I velocity structure into the H I velocity structure. Once this 
conversion is completed, the D/H ratio in the Complex C components can be estimated by min- 
imizing the residual absorption in the negative velocity wings of the H I Lyman-series lines. To 
provide the most accurate H I velocity structure determination, we placed several constraints on 
the H I model profiles constructed with the O I template, as follows: 

1. We required that the column densities of the three principal component groups (Complex C, 
IV Arch, and ISM) have total H I column densities determined from the H I 21 cm emission 
data. The H I column density of a fourth component group, the positive IVC, was calculated 
from N(0 I) assuming a solar abundance of oxygen. These column densities are summarized 
in Table 6. Within these component groups, we distributed the H I column densities in the 
proportions derived for the O I components (see Table 4). For example, in Complex C, we 
required that the -127.4 kms~^ component contain ~ 97% of the total H I column density. 
We call the column density fraction in the Complex C components xc, and set xci = 0.97 
and = (1 ~ ^Ci) = 0.03. 

2. We required that the adopted H 1 model reproduce the available 21 cm emission profiles. For 
Complex C, this required a match to the interferometric data described in §2.3. For the 
remaining components at vlsr > —100 kms~^, this required a match to the single-dish (9^7 
beam) Effelsberg data. Model parameters constrained by this criterion can be found in the 
notes for Table 6. 

3. We required the eight H I components to have velocities within a few kms~^ of those of the 
O I components. It was necessary to deviate from the O I velocities slightly for only two of the 
eight components. For the highest negative velocity component in Complex C (component 
#1), the difference was 2 kms^^; this deviation is well within the uncertainty of the O I 
model velocity but was necessary to reproduce the H I 21 cm emission profile. For the ISM 
feature at -6.5 km s"-*^ (component #6), the deviation was 4 kms~^. In this case, the precise 
central velocity is not well constrained since additional sub-components (which do not affect 
our analysis) may be present. The O I and H I velocities of all other components in the model 
were the same. 

4. With a few notable exceptions discussed below, we adopted the O I b-values for the H I 
components under the assumption that turbulent broadening dominates the line widths. In 
most cases the H I components are strongly saturated and overlapping, so a precise estimate 
of the individual component widths was not possible. However, the total widths of the 
component groups was well constrained by both the O I absorption-line data and the H I 
21 cm emission data. For the Complex C and positive IVC components, the widths could be 
constrained by the steep absorption walls of the higher order Lyman-series lines. For these 
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cases, the differences in the H I and O I b-values can be understood if there is a thermal 
contribution to the derived hne widths. 

5. We fixed the D/H ratio in all components except the two Complex C components at the local 
ISM value of D/H = 1.5 x 10~^ (see Moos et al. 2002 and references therein). The exact choice 
of D/H for the ISM and intermediate- velocity components does not influence the results of 
this study in any way, as long as the value chosen is not absurdly large. Our study of the 
H I lines does not provide any meaningful constraints on those D/H values; by adjusting the 
H I parameters slightly, values of D/H as large as those toward 7^ Vel (D/H = 2.2 x 10~^; 
Sonneborn et al. 2000) could be accommodated easily, as could values as low as D/H = 0. 
For the two Complex C components, the D/H ratio was allowed to vary under the constraint 
that D/H was equal in the two components. 

Terrestrial H I airglow is a major contaminant of the observed H I Lyman-series profiles. 
The width of the airglow emission lines is approximately 100 kms~^ since the airglow fills the 
FUSE LWRS apertures. Fortunately, the intensity of the emission is reduced substantially in the 
night-only data considered here, and the velocity of the H I airglow emission is sufHciently low 
that it does not occTir at the velocities expected for the D I absorption in Complex C (roughly 
-250 to -180 kms^^ with respect to the rest velocity of the H I lines). Furthermore, the airglow 
strength diminishes in the higher-order lines in the Lyman series. Since the airglow contamination 
occurs at velocities where the H I lines are expected to be strongly saturated, the airglow does not 
compromise the H I column density and velocity analyses. 

In Figure 8 we show the absorption produced by the H 1 model assuming (D/H)compiex c = 0- 
The wavelengths of the eight H I and D I components are marked above each spectrum, and 
a velocity scale for the H I lines is given at the top of each panel. The data are plotted as 
histogrammed lines. The model fit shown with the heavy dashed lines in Figure 8 contains H I, D I, 
and O I; no lines of other species are present in these panels. The contribution of the O I absorption 
to the model is shown by the smooth solid line. For the Ly6 and Ly9 lines, O I absorption occurs 
at wavelengths near the H I absorption but does not overlap the H I and D I lines. For the Lye 
line, O I absorption overlaps the low velocity portions of the profile but does not occur at the 
wavelengths expected for D I absorption in Complex C. The residual absorption in the negative 
velocity wings of the Ly5 and Lye profiles in Figure 8 is due to the absence of D I at Complex C 
velocities in the model illustrated. 

We have isolated the H I absorption due to different component groups in Figure 9, where 
we show the Complex C, IV Arch, and ISM/positivc IVC absorption in the H I LyO line. Note 
that the ISM and IV Arch produce no absorption at the expected velocities of the D I lines in 
Complex C. The steep absorption wings in the higher order Lyman-series lines set limits on the 
widths of components 1 and 2 (Complex C), and component 8 (the positive IVC). These widths 
and their associated uncertainties are listed in Table 6. We adopt a line width b = 9.6 it 1.2 kms~^ 
for component #1 based on considerations of the H I 21 cm emission. 
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We investigated the effects of adding to our model a broad H I feature associated with the 
hot gas traced by the O VI absorption in Complex C. Scmbach et al. (2003) find log N(0 VI) 
13.72ib0.13 with a central absorption velocity of -110 kms~^. The H I column density associated 
with this hot gas is 



^^^^)-(Z/Ze) (O/H)^ fo^' 
where fni and fovi are the ionization fractions of H I and O VI, respectively. 

For gas at T = 2.8 x 10^ K, the temperature at which O VI peaks in abundance in collisional 
ionization equilibrium, fni = 1.65 x 10~^ and fovi = 0.22 (Sutherland Sz Dopita 1993). The solar 
abundance of oxygen is (O/H)© = 4.90 x 10"^ (Allende Prieto et al. 2001). Thus, for Z/Zq = 0.1, 
we find N(H I) f« 8 x 10^^ cm^^. At T = 2.8 x 10^ K, the H I lines arising in the hot gas would 
have b ~ 68 kms~^ and would produce negligible absorption at the velocities of the D I lines in 
Complex C. (The maximum H I absorption depth would be ~ 1% of the continuum level in the Ly/3 
line, and ~ 0.01% in the higher order Lyman-series lines we are interested in here.) The assumption 
of ionization equilibrium in Eq. (2) is probably overly simplistic, and more realistic estimates would 
treat the (presently unstudied) non-equilibrium cooling of the gas. The H I estimates increase if 
the gas is not in collisional ionization equilibrium and is at lower temperatures than implied by 
the presence of O VI. However, for the H I and D I lines of interest in this study, the expected H I 
absorption depths at the velocities of the higher order D I absorption lines are still less than a few 
percent of the local continuum, even if the value of fn i / fo vi in the above equation is increased by 
a factor of 50-100. A similar conclusion is reached for the H I associated with the O VI absorption 
produced by the Galactic disk and halo along the sight line. Even though the Galactic O VI column 
density is higher than it is in Complex C [log N(0 VI) = 14.22 ± 0.05; Savage et al. (2003)], the 
velocity is lower and the impact on the D I absorption in Complex C is negligible. We conclude 
that the presence of H I associated with the hot gas traced by O VI in Complex C does not affect 
the results of our study. 

A check on the overall validity of the velocity structure along the sight line is provided by the 
H I and D I Lya and Ly/3 lines shown in Figure 10. These lines were not included explicitly in the 
derivation of the H I and D I models. The ability of the model to reproduce the damping wings 
of the Lja and Ly/3 absorption features so reliably in Figure 10 indicates that the overall column 
density estimates for the two dominant components along the sight line (Complex C and the ISM) 
are reasonable. 

A second check on our H I model comes from the 21cm emission data shown in Figure 11, 
where we have plotted a composite profile using the interferometric results for Complex C (—150 < 
■^LSR < —100 kms""*^), and the single-dish Effclsberg data for the IV Arch and ISM emission features 
at vi^sR > —100 kms~^. We converted the H I absorption model into an emission model under the 
assumption that the gas is optically thin, such that N(H I) = / lSi(v)dv = 1.823 x 10^^ / Th(v)dv 
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(Spitzer 1978). The agreement of the model and the 21cm data in Figure 11 is excellent, as 
required by the constraints imposed above. The weak H I absorption components (^^5 and ^1) at 
—29 kms~^ and +21 kms~^ do not produce significant 21cm emission. 

We note that the interferometric data for the lower velocity H I features toward PG 1259+593 
are consistent with the single-dish measurements, but are much noisier. The interferometer data 
show that there is relatively little small-scale structure in the IV Arch and the low-velocity gas. 
The interferometer observation automatically filters out the larger scale (> 20') features, leaving 
little signal in the ISM or IV Arch features and implying that arcminute-scale variations are small, 
no more than 10-20% of the total column density. This can be explained by the fact that the 
IV Arch and the low-velocity clouds along the PG 1259+593 sight line are relatively close and do 
not exhibit any bright compact cloud cores. 



7. Deuterium 

7.1. Complex C Profile Fitting Results 

Using the information derived above for O I and H I, it is possible to bracket the allowed 
range of values for N(D I) in Complex C. We considered several approaches for estimating the D I 
column density. Since most of the D I lines are weak, the strongest constraints on N(D I) come 
from the parameters assumed for the stronger Complex C feature (component #1). The weaker 
feature (component #2) does not contribute significantly to the observed D I absorption features 
unless the D/H ratio in this component is an order of magnitude or more higher than the D/H 
ratio in the stronger component. For simplicity, we assume that the D/H ratio is the same in both 
components, and focus our attention on the -127.4 kms~^ component. 

The widths of the O I and H I lines in the stronger component bracket the width expected 
for D I (i.e., b = 6.0-9.6 kms"^). If we assume that the differences in the b-values between O I 
and H I are due to a combination of turbulent and thermal broadening of the lines, then we can 
calculate the component widths assuming that b^ = + bthcrm- Solving this equation for the 

observed O I and H I b-values (see Tables 4 and 6) simultaneously leads to bturb = 5.7 + 1.1 kms~^ 
and T = 3600 + 1700 K. Using these parameters implies b 8 ± 1 kms~^ for D I. 

The model found by minimizing the residuals of the fit for the H I and D I lines has N(D I) 
= 2.35 X 10^^ cm-2 and {V) /B)^^^^^^^ ^ = 2.Q x lO'^. This result is insensitive to the width 
chosen for the weaker Complex C component, which was assumed to be 15 kms~^; values of b 
for component 2 spanning the widths implied by the O I and H I lines (b = 10 — 24 kms~^) are 
plausible and indistinguishable. Parameters for the D I components in this model can be found in 
Table 7. For the ISM and IVC lines, we have adopted the component widths and velocities found 
for H I, and we have set the D/H ratio to D/H = 1.5 x 10"^. 

Model results for the H I, D I, and O I absorption along the sight line are shown in Figures 12 
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and 13 for the FUSE SiC2 and SiCl data, respectively. The absorption features produced by the 
Lyman-series lines of LyJ-LyK arc shown as histogrammed lines. The total (H I + D I + O I) 
absorption model is shown as thick dashed curves. The smooth solid line isolates the contribution 
to the fit produced by O I absorption along the sight line. The velocity scale applies to the H I line 
shown in each panel. In this model, (D/H)compiexC = 2.6 x 10~^. The composite model fits the data 
reasonably well, with the noticeable exception of a small underestimate in the LyS and Lye lines 
near vlsr = —190 knis~^ on the H I velocity scale (~ —108 kms~^ on a D I velocity scale). The 
underestimated absorption at these velocities, which is highlighted with arrows in the top panels of 
Figures 12 and 13, is important because it impacts the quality of the fit for the higher-order lines. 
The model is driven toward higher values of D/H to fit the Ly5 and Lye lines at the expense of 
slightly over-producing the absorption in the weaker lines. Before adopting a final column density 
and error for N(D I), we consider some reasons for this underestimate. 

First, we considered the possibility that the underestimated absorption may be due to how the 
total D I column density in Complex C is distributed between the two components. In our model we 
have assumed that the column densities of O I, H I, and D I are divided among the two components 
such that xc\ = 0.97 and xc.2 = 0.03. Holding other parameters fixed, we allowed the ratios of 
N(D I) and N(H I) in the two components to vary. Values of xc^ as low as 0.94 improved the fit 
to the data. Values of xc^ < 0.94 produced noticeable discrepancies with the negative velocity 
absorption walls of the higher order H I lines, in the sense that too much absorption was predicted. 
Lowering the b-value of component 2 to alleviate this concern reintroduced the discrepancy in the 
LyS and Lye lines. The changes considered affected both N(D I) and N(H I) by the same percentage 
since the D/H ratio was assumed constant in Complex C. Thus, while the fit could be improved 
slightly, such changes had no effect on (L)/H)p^^pj^^ ^. 

We also considered whether a difference in the D/H ratio in the two Complex C components 
might lead to an underestimate of the absorption at lower velocities. To account for the extra 
absorption by increasing the D/H ratio in the weaker feature (component #2 at -112.5 kms~^), 
requires a value of D/H that is at least a factor of 8-10 times greater than D/H in component 1. 
Such a discrepancy is inconsistent with Big Bang nucleosynthesis estimates of the primordial value 
of D/H, (D/H)p (2.6 - 3.0) x 10"^ (NoUctt & Buries 2000; Buries et al. 2001; Spergel et al. 
2003), and the expected variations in D/H caused by deuterium astration. 

Next, we explored the possibility that the residual absorption could be due to an underestimate 
in the width of the D I and H I lines in the main Complex C component at -127.4 kms~^. To fit 
the absorption well would require a width bni ~ 14 kms~^, which is inconsistent with the width 
of the H I 21cm emission shown in Figure 11. It is possible for the H I absorption width to be 
larger than that found for the 21 cm emission if the wings of the absorption consist primarily of 
low column density gas that is below the detection limit of the radio observations. In this case, the 
D I line width could also be increased, resulting in a slightly better fit to the data. The resulting 
D/H value is similar to that found for the single-component curve of growth analysis with boi ~ 14 
described in the next section. 
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The underestimated absorption may be due to a small amount of fixed-pattern noise in the 
data for one or both of the SiC channels. Note that the underestimate is slightly different for the 
lines shown in Figure 12 for the SiC2 channel and those shown in Figure 13 for the SiCl channel. If 
the absorption is actually fixed-pattern noise, the fit may change the derived value of N(D 1) since 
the D I profile fit is sensitive to absorption at this velocity. 

We considered whether adding an additional component to the absorption model to account 
for the residual absorption affected the fit results. This feature could be either a D I component 
corresponding to a stronger H I feature near -108 kms~^, or a weak interloping H I component 
near -190 kms~^. To account for the absorption with additional D I requires an H I component 
near -108 kms'^ with N(H I) « (2 - 5) x 10^^ cm^^ if D/H = (1.5 - 3.0) x 10"^. While this 
additional H I component would remain hidden in the H I Lyman-series lines, it would lead to a 
serious confiict with both the observed H I 21 cm emission and the O I absorption at this velocity. 

An H I feature at -190 kms^^ with N(H I) 3 x lO^'' cm'^ and bm « 10 - 15 kms'^ 
can account for the absorption and remove the noted discrepancy. The weak high-velocity H I 
interloper would be undetectable through either its 21cm emission or corresponding metal- line 
absorption (see §6.1). The associated D I absorption would not be detectable at higher velocities. 
In considering this possibility, we note that the Lya and Ly/3 lines do not help constrain whether 
the absorption could be H I because in both cases the features occur in the strong cores of the 
overall sight line absorption, hyj is also not restrictive since strong O I absorption blends with 
D 1 and H 1 absorption at this velocity. If we add a weak H I feature with these properties to our 
model and recalculate the fit, we find that the preferred value of N(D I) in Complex C drops from 
2.35 X 10^5 cm-2 to 1.62 x lO^^ cm-^, and (D/H)^^^^,^^ ^ drops from 2.6 x 10"^ to 1.8 x 10"^ 
We reproduce Figures 12 and 13 in Figures 14 and 15 with the additional weak H I feature added 
to the model (i.e., D/H = 1.8 x 10~^). The greater the contribution of this feature, the smaller 
the inferred D/H ratio in Complex C becomes. Many additional parameters could be tuned to 
improve the fits shown in Figures 14 and 15, but the degeneracies in the possible solutions preclude 
a definitive statement about the exact magnitude of the absorption caused by this weak feature. 
The discrepancies in the fit caused by the presence or absence of this feature lead to an uncertainty 
of ±0.4 X 10^^ cm~^ in the D I column density. 

It is important to note that H I 21 cm emission with N(H I) ~ 10^^ cm~^ is seen near -190 
kms~^ along a high-velocity ridge running through Complex C (see Figure 6 in Wakker 2001 and 
Figure 13 in Sembach et al. 2003). PC 1259+593 lies off of this ridge by several degrees, and even 
though there are no metals detected in our absorption-line spectra of PC 1259+593, there is still a 
reasonable possibility that there may be low column density H I near -190 kms~^ in this direction. 
In our final estimation of the deuterium abundance in Complex C, we will account for the possible 
range of values of (D/H)compiex c resulting from the inclusion or exclusion of a weak H I feature. 
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7.2. Curve of Growth Results for Complex C 

As an additional consistency check on our profile fitting results for the D I column density in 
Complex C, we constructed curves of growth for the Complex C absorption using the equivalent 
widths of the D I lines measured between —160 < vlsr < —110 kms~^. To make the D I equivalent 
width measurements listed in Table 8, we removed the H I and O I signatures in the data by dividing 
the FUSE spectra by the adopted O I and H I models with the parameters listed in Tables 4 and 6. 
Next, we measured the D I absorption in the 50 kms~^ interval in each channel between -160 
and -110 kms~^ (-242 to -192 kms^^ on the H I velocity scale shown in Figures 12-15). This 
velocity range, which is indicated by the horizontal bar above the Lye spectra in Figures 12 and 
13, includes the Complex C D I absorption recovered by our model result but does not include the 
small amount of residual absorption at lower velocities discussed in the preceding section. 

The results of this process are shown in Figure 16, where we plot the data along with both a 
single-component curve of growth and the double-component curve of growth based upon the D I 
fit parameters listed in Table 7. The best-fit single-component curve of growth yields log N(D I) = 
15.22ibQj2 (lo") and b = 14.4ib4^5^ kms^"*^ (l*^)- This column density is less than that implied by 
the profile fitting process (at least for the case where there is no H I interloper) and the optimum 
b-value is higher than the H I b-value for the dominant component. However, it is still within 
the possible ranges allowed by the uncertainties in the profile fits. The single-component COG 
b-value is determined in large part by the measurements for the hy5 line in the two channels. Since 
the measurements for this line differ between the two channels by an amount similar to the la 
measurement uncertainties, the resulting b-value is particularly uncertain. For comparison, the 
dashed curves in the top panel of Figure 16 show the single-component COGs for b-values of 10 
and 20 kms~^, assuming the same value of N(D I). 

If we take the additional residual absorption redward of the Complex C absorption into account 
and assume that it is all due to D I, the column density from the single-component curve of growth 
method becomes log N(D I) = IB.SSztQ-Jf and b = 26±Qg kms~^. We conclude that the quality of 
the fit docs not improve, which is consistent with our contention in the previous section that this 
residual absorption is unlikely to be due to D I. 

The double-component COG shown in the bottom panel of Figure 16 passes through the SiCl 

Lj^^ measurement but underestimates the SiC2 LyS measurement slightly, as expected from the 
profile fit results shown in Figure 12 and 13. The SiC2 equivalent width measurement includes 
some residual absorption located on the sides of the Complex C absorption that is not fit well by 
the model. This discrepancy between the LyS measurements for the two channels may be due to 
fixed-pattern noise in the data (see discussion of errors below). 

The curve-of-growth results illustrate an important point in the consideration of the Complex C 
D I column density. The predicted strengths of the lines are sensitive to the value of bni chosen 
for the stronger component because the LyS line does not lie on the linear portion of the curve 
of growth. Unfortunately, omitting the LyS measurements is a poor alternative to including them 
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since the remaining detections are weak and have relatively large errors. Thus, it is important to 
consider the systematic error associated with a range of possible b-values. 

7.3. Sources of Systematic Errors 

In addition to our concerns about the ability of the profile fitting model to reproduce the small 
amount of residual absorption in the H I + D I absorption profiles, there are other sources of error 
that must be considered in making an estimation of N(D I) in Complex C. The dominant systematic 
errors include: 

1. Inability to distinguish between small column density differences due to the unresolved nature 
of the strongest Complex C component. This is probably the most troublesome problem since 
it necessarily guarantees that the minimum in N-b parameter space is shallow. This occurs 
because the main D I component in Complex C has a width (FWHM ~ 15 kms~^) that is 
less than the width of the FUSE LSF (FWHM 20 - 25 kms-^). Using the O I and H I 
widths as guides, we estimate that the error arising from the uncertainty in b when holding 
other parameters fixed is 0.4 x 10^^ cm~^. 

2. The assumed velocity structure for the sight line, particularly the Complex C and intermediate 
velocity components. The error associated with the assumed velocity structure is difficult to 
quantify. One important example, that of the possible existence of very high velocity H I 
absorption, was provided above and must be considered in the error budget (see point 4 
below). It is also possible that the velocity structure for the stronger Complex C component 
consists of several narrow absorption features closely spaced in velocity. High velocity clouds 
typically consist primarily of warm gas and do not often exhibit evidence for cold components 
(Wakker & van Woerden 1997). Most HVCs show no evidence of molecular gas, including 
Complex C (Richter et al. 2001a). In one region of the Magellanic Stream containing dust in 
which H2 has been detected in absorption by FUSE (Sembach et al. 2001a), the cloud structure 
consists of a compact cloud core that is surrounded by warmer gas. The high-resolution STIS 
data for Complex C do not show any evidence of C I absorption that might be associated 
with cold components in the PG 1259+593 direction. Therefore, if narrow components do 
exist, there is probably no single cold component that dominates the absorption. In this case, 
the single-component approximation usually results in column densities accurate to better 
than 20% (Jenkins 1986; Savage & Sembach 1991), or about 0.5 x 10^^ cm"^ ^j^^ j 
Complex C absorption toward PG 1259+593. An additional source of uncertainty related to 
the Complex C velocity structure results from the component width and strength trade-offs 
that can be made between the strong and weak components (components #1 and #2). We 
find that these typically affect the total column density at a level less than 0.2 x 10^^ cm~^. 

3. Fixed-pattern noise that may he present in these FUSE data. Small differences are visible 
in the SiC2 and SiCl data for some of the profiles considered in this work. For example, 
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the SiCl data for the Ly5 line recover fully to the continuum blueward of the Complex C 
absorption, while the SiC2 data do not (see the top panels of Figures 12 and 13). There are 
also slight discrepancies in the wings of the two hj5 lines, which lead to the differences in D I 
equivalent widths listed in Table 8. If we consider the difference in column densities resulting 
from the omission of the SiC data for one LjS line or the other as a measure of the error 
caused by fixed-pattern noise in these lines, we find an uncertainty of 0.5 x 10^^ cm~^ in 
the D I column density. 

4. The possibility that the absorption we identify as D I in Complex C might be high velocity H I. 
We considered the case of the weak residual absorption above. However, our analysis of the 
metal lines and the inferred limits for the amount of H I at high velocities cannot formally 
exclude the possibility that some of the other absorption is high velocity H I. If all of the 
observed absorption were H I and not D I, the D/H ratio in Complex C would be strange 
indeed! We acknowledge the possibility that some of the absorption assumed to be due to 
D I in Complex C could be due to weak, high-velocity H I features. However, without further 
evidence of this, we assume that the absorption observed is caused by D I, with the caveat 
that the D I column density decreases by ~ 0.8 x 10^^ cm~^ if a weak H I feature is included 
to explain the residual absorption near the D I absorption. 



7.4. Complex C Column Density 

We adopt a D I column density for Complex C that is an average of the values of N(D I) for 
the case where the residual absorption is left untreated in our model (no H I interloper, N(D I) 
= 2.4 X 10^^ cm^^, as shown in Figures 12 and 13) and the case where the residual absorption is 
assumed to be due to an H I interloper (N(D I) = 1.6 x 10^^ cm~^, as shown in Figures 14 and 
15). To derive an error, we add in quadrature the uncertainties associated with this averaging and 
the systematic errors described in §7.3. We treat the individual uncertainties as 2(7 estimates and 
combine them in quadrature because the errors do not all add with the same sign. We find 

N(D I) = (2.0 ± 0.6) X 10^5 cm-2 (68% confidence) 
N(D I) = (2.0 ± 0.9) X 10^5 cm-2 (95% confidence), 

where we have scaled the 2a uncertainty by the confidence interval to produce the la uncertainty. 
The upper 2a limit produces a noticeably unsatisfactory fit to the data. 



8. Discussion 



We summarize the Complex C column densities of H I, D I, and O I in Table 9, together 
with the la (68% confidence interval) and 2a (95% confidence interval) error estimates based on 
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the descriptions of errors given above. Systematic errors dominate the uncertainties for all three 
species. We also list the D/H, D/0, and 0/H cohimn density ratios, along with the propagated 
errors from the column density determinations. A summary of the light element abundance ratios 
for Complex C and other environments is given in Table 10, along with references for the ratios. 
The environments considered in this table include the local ISM, the ISM of the Galactic disk at 
distances of 200-1000 pc from the Sun, and high-redshift absorption line systems (both Lyman- 
limit systems and damped Lya systems). We list both D/H and D/0 when estimates of [0/H] 
are available. These comparisons are not meant to be exhaustive summaries of all the available 
data in the literature; rather they are synopses of several studies given to provide insight into 
how the Complex C values compare to typical values found elsewhere. The data in Table 10 are 
shown graphically in Figure 17, where we plot D/H as a function of metallicity determined through 
measurements of [0/H], or [Si/H] if [0/H] is not available. For the Lyman-limit systems where 
[Si/H] is used, the metallicity may be somewhat uncertain because large ionization corrections 
are necessary to convert estimates of N(Si II)/N(H I) into estimates of N(Si)/N(H). Various other 
metallicity indicators are available, but [0/H] tends to be the best for the reasons described in 
§5 (see also the detailed discussion presented by Timmes et al. 1997). Our estimate of D/H in 
Complex C rules out values of D/H greater than 3.3 x 10~^ or less than 1.1 x 10~^ with a reasonable 
degree of confidence (~ 95%). Higher values over-produce the amount of D I absorption in the 
higher order Lyman series lines, and lower values under-produce the amount of absorption expected 
in the hyS line. Both limits have interesting implications. 

Most models of Galactic chemical evolution predict that astration of deuterium should result in 
a factor of ~ 1.5 — 3 decline in the cosmic abundance of deuterium from the Big Bang to the present 
time (Edmunds 1994; Tosi et al. 1998; Chiappini et al. 2002), though higher levels of astration may 
be possible if special conditions, such as a Galactic wind, are invoked to counter over-production 
of elements heavier than helium (see, e.g., Scully et al. 1997). For a simple closed-box model with 
the assumption of instantaneous recycling, a deuterium astration efficiency of 60%, and standard 
oxygen yields, Pagel (1997) finds that the there should be only a few percent reduction of deuterium 
relative to its primordial abundance if [0/H] < — 1, as found for Complex C. In Figure 17, we show 
the results for a simple chemical evolution model from Fields et al. (2001) having (D/H)p and 
little destruction of deuterium at low metallicities from Population III stars. More complicated 
chemical evolution models requiring bimodal episodes of star formation to explain the scatter in 
D/H values measured at high redshift can be found in that work. Our 2a upper limit of D/H 
< 3.3 X 10~^ suggests that the value of D/H produced by Big Bang nucleosynthesis cannot be 
too much larger than this upper limit unless the assumptions about the destruction of deuterium 
and production of oxygen are invalid. Similar arguments have been made to infer a primordial 
value of (D/H)p < 3 x 10~^ from determinations of D/H in high-redshift, low-metallicity systems 
(Tytler et al. 2000; O'Meara et al. 2001; Kirkman et al. 2003). Values of (D/H)p < 3 x lO^^ 
imply values of the baryon density, Q^h^, and the baryon-to-photon ratio, 77, in concordance with 
estimates derived from cosmic microwave background (CMB) measurements. Using the standard 
Big Bang nucleosynthesis predictions for the abundance of deuterium (Buries et al. 2001), a value 
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of (D/H)p < 3 X 10-^ implies n^h'^ > 0.02 and > 5.6 x 10-^°. The CMB estimates cluster near 
values of n^h'^ w 0.022 - 0.024 (Netterfield et al. 2002; Pryke et al. 2002; Spcrgcl et al. 2003), 
and imply (D/H)p = (2.62±[j:^^) x 10"^ (Spergcl et al. 2003; see also Steigman 2003). The region 
between the dashed horizontal lines in Figure 17 is the ±2c7 range of (D/H)p allowed by the recent 
Wilkinson Microwave Anisotropy Probe (WMAP) observations. 

With the PG 1259+593 data, we can confidently rule out very high D/H values for Complex C, 
like those claimed for some intergalactic systems [e.g., D/H = (2.0 ± 0.5) x 10~^ at z pa 0.7 toward 
PG 1718+4807; Webb et al. 1997]. This is important because the Complex C measurement is the 
only other measurement of D/H outside the Milky Way at 2; < 1.0. Such high values are also 
inconsistent with the CMB measurements, and recently the value for the z = 0.7 absorber toward 
PG 1718+4807 has been challenged on other grounds (e.g., Kirkman et al. 2001 - but see Crighton 
et al. 2003). 

On the low end of the D/H range, we cannot exclude the possibility that D/H in Complex C is 
similar to D/H in the local ISM (D/H ~ 1.5 x 10"^ Moos et al. 2002; Linsky 2003). However, we 
can rule out D/H values as low as those found for some Galactic disk sight lines extending beyond 
the local ISM (i.e., at c/ > 100 pc). Values of D/H less than 1.0 x 10~^ have been determined for 
two extended sight lines in the Galactic disk using FUSE data. Hoopes et al. (2003) find D/H = 
(0.85±[]J^) X 10"^ (2cr) toward HD 195965 (d ~ 800 pc), and D/H = (0.78±[j:i) x 10"^ (2cr) toward 
HD 191877 {d ~ 2200 pc). Using IMAPS data, Jenkins et al. (1999) find D/H = (0.74±[]:i|]) x 10^^ 
(90% confidence) toward 6 Ori A (d ~ 500 pc); a similar result was found by Laurent, Vidal-Madjar, 
&; York (1979) using Copernicus data. The low values found for these sight lines indicate that more 
astration of deuterium may have occurred in these regions than in the local ISM or in Complex C. 
The local star formation histories of the regions explored toward S Ori A and HD 195965 are roughly 
consistent with this idea (see Hoopes et al. 2003 and Jenkins et al. 1999), but the number of sight 
lines is still too small to draw general conclusions. Eventually, with a large enough number of sight 
lines, it may be possible to determine whether refinements to chemical evolution models are needed. 

The value of D/0 = 0.28 ± 0.12 in Complex C is very different from the D/0 ratio in the 
local ISM or the disk of the Milky Way. Hebrard & Moos (2003) find D/0 = 0.038 ± 0.002 for 
white dwarf and subdwarf sight lines within ~ 150 pc of the Sun. Most of this difference can 
be accounted for by the different metallicities of the local ISM and Complex C. Note that even 
without the H I measurement, the D/0 ratio strongly suggests that the metallicity in Complex C 
is lower than solar. The Complex C result is more similar to the value of D/0 = 0.37 ± 0.03 found 
for the z = 3.025 absorber toward Q 0347-3819 (Lcvshakov ct al. 2002) than it is to the value 
for the Galactic disk. However, values of D/0 for two other high-redshift systems are an order of 
magnitude or more larger than the Complex C value. We note that determinations of N(0 1) in the 
high-redshift systems are challenging since the weaker O I lines occur within the H I Lya forest. 

As estimates of D/H in the high-redshift systems have fallen in recent years, a few local 
measurements with high values of D/H have received more attention. For example, the value 
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of D/H = (2.18±g|?) X 10"^ for the 7^ Vel sight hne (90% confidence; Sonneborn et al. 2000) 
is indistinguishable from the vahies for the high-redshift (z > 2) systems toward HS 0105+1619 
and Q 1243+3407, and pcrliaps liigher than the value for Q 2206-199 (see Table 10 and references 
therein). Preliminary FUSE results for the IX Vel sight line, which lies ~ 1° from 7^ Vel, indicate 
that the value of D/H for the first ~ 100 pc of the 7^ Vel sight line is similar to that in the local 
ISM (Blair et al. 2003). If this result is confirmed by upcoming FUSE observations, then the D/H 
ratio in the Vela region must be even higher than the sight line average. Note, however, that the 
Pup sight line, which is in the same region of the sky and extends further than the 7^ Vel sight 
line has a D/H ratio similar to that of the local ISM (see Table 17). One possible explanation for 
a high value in the ISM toward 7^ Vel is infall of "primordial" (metal-poor, D-rich) gas. Clouds as 
large as Complex C could contribute significant amounts of deuterium to regions of the Galactic 
disk; a mass of M(H I) > (1.2-3.0) x 10^ Mq and D/H = 2.2 x 10"^ imply that Complex C contains 
> (26 — 66) M.Q of deuterium. The effects of infalling low-metallicity gas have been considered by 
various authors in the context of deuterium astration and potential solutions to the well-known 
"G-dwarf" problem (Tosi 1988; Mattcucci & Francois 1989; Tosi et al. 1998; Wakkcr et al. 1999). 
However, a serious problem with this interpretation for the 7^ Vel region is that the gas along the 
sight line does not show a corresponding reduction in the ratios of metals to hydrogen as would 
be expected if infall were the explanation for the high D/H ratio. Furthermore, the constancy 
of the oxygen and krypton abundances within several hundred parsecs of the Sun (Meyer, Jura, 
& Cardelli 1998; Cartledge et al. 2001) supports the idea that the nearby gas is reasonably well 
mixed and has similar heavy element abundances. Any explanation for the high D/H ratio in the 
Vela region that incorporates infall must also include processing of the gas to explain the relatively 
"normal" metal abundances in the region. 

Other possibilities also exist for explaining the variations of D/H in the Milky Way. For exam- 
ple, D could be depleted onto dust grains (Jura 1982). From simple thermodynamic considerations, 
Draine (2003) has shown that differences in the binding energies of D and H on the peripheries 
of polycyclic aromatic hydrocarbons (PAHs) could produce extreme enrichment of deuterium in 
bound form, similar to what has been found for some simple interstellar molecules. He argues 
that it is plausible that roughly 20% of the H atoms bound to PAHs are replaced by D. If this is 
correct, he estimates that there are sufficient PAH binding sites in a volume of normal interstellar 
material to account for a decrease of free D atoms by one part in 10^ of H (i.e., roughly the size of 
the reduction observed along some lines of sight). In some places, the atomic medium can revert 
to the intrinsic D/H when the PAHs are destroyed by the recent passage of a shock, which could 
explain why D/H seems to vary. Reinforcing the concept that variable D abTindances might be 
explained by Draine's proposal is the recent detection of the CD bond stretch features at 4.4 and 
4.67 iim toward the Orion Bar reported by Peeters (2002). She found D/H = 0.17 + 0.07 in bound 
form. This detection, although not very far above the noise, indicates that Draine's estimate of 
a 20% replacement of H by D atoms may be approximately correct. For any gas with an overall 
metallicity as low as that of Complex C, the concentration of PAHs would be so low that there 
would be no measurable reduction in the concentration of free deuterium atoms, which is consis- 
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tent with our finding for the sight line toward PG 1259+593. Finally, Mullan & Linsky (1998) have 
proposed that deuterium production in stellar flares may cause some variations in D/H locally, but 
Prodanovic & Fields (2003) have recently reconsidered this issue and find that flare production of 
deuterium is not a significant source of D I on galactic scales. They concur with Mullan & Linsky 
that local variations by stellar flare production cannot be ruled out. Dust depletion and stellar flare 
production of deuterium are unlikely to alter the present value of D /H in Complex C substantially. 
Complex C contains no known stars, and analyses of its metal content indicate that it contains 
little dust. 

The metallicity of ZjZ^ 0.17 we derive for Complex C through the ratio of O I and H I is 

very similar to the recent determination of the metallicity in Complex C toward 3C 351 by Tripp 
et al. (2003). They find [0/H] = — 0.76±Q2f. Both metallicity estimates are slightly higher than 
the value of Z/Zq ~ 0.1 inferred for the Mrk 290 line of sight through Complex C based on 
measurements of S II and H I (Wakker et al. 1999). The metallicity estimates may depend on 
environment and may change slightly as a function of position within Complex C; Collins et al. 
2003 find a range from 0.1 to 0.25 solar, but some of these estimates are complicated by uncertain 
ionization corrections for S II. One possible explanation for the slightly varying abundances in 
Complex C might be that the gas is interacting with the Galactic corona or high Galactic halo, as 
indicated by measurements of highly ionized gas associated with the complex (Sembach et al. 2003; 
Fox et al. 2003). Mixing of the diffuse regions of Complex C with the more metal abundant gas of 
the Galactic thick disk and halo would tend to increase the observed abundances in these regions, 
while the abundances in denser clumps of gas, as seen toward PG 1259+593, should reflect more 
closely the original (pre-inter action) abundances in the cloud. 

It is interesting that our best-fit value of D/H found for Complex C is similar to that found 

for high-redshift damped Lya systems (see Pettini & Bowen 2001 and Kirkman et al. 2003). The 
H I column density of N(H I) = 9.0 x 10^^ cm^^ is sufficient to produce damping wings on the 
H I Lya and Ly/3 lines. If observed in isolation from the Milky Way by an outside observer. 
Complex C would produce H I and D I absorption similar to that seen in these higher redshift 
systems. Complex C also has a low nitrogen abundance inferred from the ratio of N I to H I, [N/H] 
« -1.9 toward PG 1259+593 (Richter et al. 2001b; Collins et al. 2003), about a factor of 10 lower 
than [0/H]. (Initial work on the amount of N II present indicates that photoionization of N I in 
the H I regions is not able to account for the deficit in the derived value of [N/H].) Commensurate 
values of [N/H] have been found for other Complex C sight lines (Collins et al. 2003; Tripp et 
al. 2003). Some damped Lya systems show similar deficits of nitrogen relative to oxygen, which 
several authors have explained as a nucleosynthetic effect resulting from delayed release of nitrogen 
between massive and intermediate mass stars or a top-heavy or truncated initial mass function 
(Prochaska et al. 2002; Molaro 2003). At the metallicity of Complex C, the low N/0 ratio observed 
is probably more dependent upon the secondary nature of nitrogen production than these other 
processes (see Figure 6 in Pettini et al. 2002). The metallicity, abundance pattern, D/H ratio, 
and lower distance limit all indicate that Complex C is an external system (or the remains of an 



-35- 



external system) falling into with the Milky Way rather than gas ejected from the Galactic disk. 

9. Future Progress 

Given the complex velocity structure of Galactic lines of sight, it may be difficult to find 
other HVCs for which an analysis similar to the one presented here can be performed. A similar 
conclusion was reached by Hoopes ct al. (2003) for Galactic disk ISM along sight lines extending 
more than ~ 1 kpc. We have inspected the FUSE spectra of several other Complex C sight lines, 
but none offer comparable prospects for a second determination of the D/H and D/0 ratios in 
Complex C. Some brief comments on these sight lines are contained in Table 11. We are currently 
inspecting data for several other HVCs. 

Reducing the error bars on the Complex C D/H determination toward PG 1259+593 is desir- 
able. A reduction in the error bars from a la uncertainty of ~ 30% to ~ 10% would allow stronger 
statements to be made about the abundance of deuterium in this environment and the evolution 
of deuterium as a function of metallicity. It would be difficult to improve upon the current values 
for PG 1259+593 with better FUSE data. A large amount of additional observing time (~ 2 x 10^ 
seconds) would be needed to reduce the statistical errors by a factor of ^ 2. This would improve the 
continuum estimations and the O I velocity model, but it is unclear whether this would translate 
into a comparable reduction in the systematic errors in the H I and D I models. 

A substantial improvement in the accuracy of the D /H determination in Complex C could be 
made with higher spectral resolution far-ultraviolet absorption-line data. Even a factor of 1.5-2.0 
increase in resolution to i? ~ 30, 000 — 45, 000 (similar to that available with cchcllc spectro- 
graphs on large ground-based telescopes or with HST/STIS at ultraviolet wavelengths > 1150 A) 
would significantly enhance the discrimination between Complex GDI and H 1 absorption in the 
Lyman-series transitions toward PG 1259+593. This is currently beyond the capabilities of exist- 
ing instrumentation, but perhaps someday it will be possible to observe PG 1259+593 at higher 
spectral resolution in the far-ultraviolet spectral region with a large space telescope equipped with 
an efficient high-resolution spectrograph. If and when that day comes, we recommend a re-analysis 
of the PG 1259+593 sight line velocity structure and a redetermination of the range of D/H values 
allowed by the data. 

10. Summary 

We report the detection of D I Lyman-series absorption in high velocity cloud Complex C, a 
low-metallicity gas cloud falling onto the Milky Way. This is the first detection of atomic deuterium 
in the local universe in a location other than the nearby regions of the Galactic disk. The primary 
results of our study are as follows: 
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1. We combine data from the Far Ultraviolet Spectroscopic Explorer, Hubble Space Telescope, 
and Westerbork Synthesis Radio Telescope to derive the column densities of H I, D I, and O I 
in the direction of PG 1259+593 (/ = 120?56, b = 58?05). 

2. We use the FUSE and HST/STIS absorption-line data to construct a velocity model for the 
sight line based on the numerous O I absorption-lines detected. We identify 8 absorption-line 
components common to all three species (H I, D I, and O I) along the sight line. Two of 
these components, at vlsr ~ —128 and -112 kms~^, are identified with Complex C. The 
remaining lower-velocity components are identified with the Intermediate Velocity Arch, the 
Milky Way ISM, and a positive intermediate-velocity cloud in this direction. 

3. Our WSRT interferometer map of the H I 21 cm emission toward PG 1259+593 indicates that 
the sight line passes through a compact concentration of neutral gas in Complex C. We use 
the WSRT data having a resolution of ~ 1' together with single-dish data from the Effelsberg 
100-mctcr radio telescope to estimate the H I column density in Complex C and to constrain 
the velocity extents of the Complex CHI Lyman-series absorption components observed by 
FUSE. 

4. We use a combination of analysis techniques to estimate the column densities of H I, D I, and 
O I in the absorbing regions along the sight line. We combine line-profile fitting analyses and 
curve-of-growth techniques with minimization codes to estimate the most likely value and 
uncertainties on the column density of each species. The use of multiple analysis methods 
provides important consistency checks on the results of the individual methods. The good 
agreement found from independent analyses by members of our team confirm that the formal 
errors are reliable representations of the uncertainties involved in the analyses. 

5. For each species, we describe the uncertainties involved in the estimation of the column 
density. Systematic errors dominate the uncertainties, especially for D I. In estimating the 
uncertainties on D/H, we have included the possibility that some of the absorption in the 
vicinity of the D I lines might be caused by a weak, high-velocity H I feature near -190 
kms"-*^. 

6. For the Complex C absorption, we find N(H I) = (9.0+1.0) x 10^^ cm-^, N(D I) = (2.0+0.6) x 
10^^ cm~2, and N(0 I) = (7.2+2.1) x 10^^ cm"^. The uncertainties quoted are la errors (68% 
confidence intervals); two sigma error estimates can be found in Table 9. The corresponding 
light element abundance ratios are D/H = (2.2 + 0.7) x 10"^, 0/H = (8.0 + 2.5) x 10-^ and 
D/0 = 0.28 + 0.12. 

7. The metallicity of the Complex C gas in this direction is approximately 1/6 solar, as inferred 
from the oxygen abundance of [0/H] = -0.79+[|:{^ (la). This metallicity is in excellent 
agreement with a recent high-quality determination of [0/H] in Complex C toward 3C351, 
which is located ~ 30° from PG 1259+593. Two other Complex C sight lines analyzed by 
Collins et al. (2003) yield similar metallicities (Mrk 279: [0/H] = -0.71+^ ^; Mrk 817: [0/H 
= -0.59+g:?f). 
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8. While we cannot rule out a value of D/H for Complex C similar to that found for the local 
ISM (i.e., D/H ~ 1.5 x 10~^), wc can reasonably exclude values as low as those determined 
recently for extended sight lines in the Galactic disk (D/H < 1 x 10"^), provided that the 
absorption observed is really D I and not very high velocity H I. This conclusion lends support 
to the hypothesis that Complex C is a manifestation of gas originating outside the Milky Way, 
rather than material recirculated between the Galactic disk and halo. 

9. A value of D/H ~ 2.2 x 10~^ for Complex C is consistent with the primordial abundance 
of deuterium inferred from recent WMAP observations of the cosmic microwave background, 
combined with simple chemical evolution models that predict the amount of deuterium astra- 
tion with time. Reducing the errors on the D/H estimate in Complex C toward PG 1259+593 
will be difficult, but may eventually be possible with a large space telescope equipped with a 
high-resolution ultraviolet spectrograph. 

10. The number of sight lines for which D I can be observed in Lyman-scries absorption beyond 
a few hundred parsecs of the Sun appears to be small because of H I velocity crowding. We 
considered the data for several other sight lines through Complex C and found that they are 
unsuitable for a reliable D I column density determination. We continue to examine FUSE 
data for D I absorption in other HVCs. 
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Fig. 1. — See caption on next page. 
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Fig. 1.— Fully reduced FUSE spectra of PG 1259+593 between 915 and 955 A. The two night-only 

datasets shown were obtained with the SiCl and SiC2 channels. The spectra are binned by 8 pixels 
(slightly less than a resolution element) for illustration purposes only. The locations of prominent 
H I, D I, N I, and O I lines are indicated at the top of the top panel. Four components corresponding 
to absorption groups at -128 kms~^ (Complex C), -54 kms~^ (IV Arch), -5 kms~^ (ISM), and 
+69 kms"^ (Positive IVC) are shown for each transition. The H I Lyman-series lines are identified 
in the lower panel. The heavy dashed lines indicate the continua adopted in our analysis of these 
spectra. The dotted line at the bottom of each panel marks the zero flux level. The wavelengths 
of terrestrial H I airglow features in these night-only observations are denoted by crossed circles. 
Note that the strong cores of the H I lines reach zero residual intensity at wavelengths that are 
unaffected by these airglow features. The source of the slight rise in the SiCl spectrum below 917 A 
is probably residual airglow at these wavelengths. 
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Fig. 2.— 



See caption on next page. 
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Fig. 2. — Top: Brightness temperature versus LSR velocity for the H I 21 cm emission observed 
toward PG 1259+593 with the Effelsberg 100-meter telescope and NRAO 140-foot telescope. The 
Effelsberg data have a beam size of 9^7 and a velocity resolution of 1.8 kms^-*^ after Hanning 
smoothing. The NRAO data have a beam size of 21' and a velocity resolution of 1 kms^^. Bottom: 
Ultraviolet absorption-line profiles of various interstellar species observed with HST/STIS. The data 
have a resolution of ~ 7 kms~^ (FWHM). The velocities of the 8 absorption components along the 
sight line are indicated above the top spectrum and are labeled C (Complex C), G (Galactic ISM), 
or I (intermediate velocity ISM). Note the multiple component nature of the Complex C absorption 
between -150 and -100 kms""*^. The weak intermediate-velocity feature at ~ -|-69 kms^^ is visible 
in the strong O I and C II lines shown; it is also visible in the H I lines shown in Figure 1. The 
Al II line falls near the edge of a STIS echelle order and extends only to -|-60 kms~^. The broad 
absorption feature between -1-100 and -1-200 kms~^ next to Si II A1526.707 is an intergalactic Lya 
absorber at redshift z = 0.2564. 
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Fig. 3. — Filtered H I column density of the WSRT data, integrated over the velocity range —148 
to —109 kms~^, which includes core CIII. There is a slowly varying offset of (3.0 — 6.0) x 10^^ 
cm~^ across the field that is not included in this estimate because it has been filtered out by the 
interferometer. At the position of PG 1259+593, marked by the star at a=13''01'"13^ (5=59°02'06", 
this offset is 6 X 10^^ cm~^. The circles centered on PG 1259+593 show the FWHM of the Effelsberg 
and Westerbork primary beams. 
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Fig. 4.— 



See caption on next page. 
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Fig. 4.— A comparison of O I lines toward PG 1259+593 observed by FUSE and HST/STIS to our 

model for the O I velocity distribution along the sight line. The data are shown as histogrammed 
lines. For A < 1000 A, the FUSE SiC2 data are shown as solid lines, and the SiCl data arc shown as 
dashed lines. The O 1 A1039.230 line is from the FUSE LiFl channel. In all cases, the model results 
are shown as smooth, heavy lines. Vertical tick marks above the STIS O I A1302.168 spectrum 
denote the velocities of the 8 components in the model. The deep broad absorption features near 
some of the O I lines are interstellar H I Lyman-series lines. The cross-hatched emission feature in 
the O I A1302.168 spectrum is terrestrial O I airglow entering the STIS slit during the observation. 
The FUSE profiles shown do not contain such features since we show only orbital night data for 
which the terrestrial O 1 emission is negligible at these wavelengths. 
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Fig. 5. — Curve of growth results for the O I absorption in Complex C between t'LSR = —180 kms^^ 
and vlsr = —105 kms~^. The x-axis has units of logarithmic inverse centimeters. Filled circles are 
measurements based on data from the FUSE SiC2 channel. Open circles are measurements based 
on data from the FUSE SiCl channel. The points for the two channels have been offset slightly 
(plus or minus 0.02 dcx) from their nominal values in the horizontal direction for clarity. The filled 
diamond indicates the HST/STIS E140M measurement for the O I A1302.168 line. Error bars are la 
estimates. Top: Single-component COG fit with corresponding profile fit to the O I A1302.168 line 
(inset). The two arrows in the inset panel indicate the velocities (in kms~^) of the two Complex C 
components. Bottom: Double-component COG fit based on the fitting parameters used to construct 
the profiles shown in Figure 4 with corresponding profile fit to the O I A1302.168 line (inset). The 
two arrows in the inset panel indicate the velocities of the two Complex C components. 
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Fig. 6. — Double-component COG fit to the O I lines based on the fitting parameters used to 
construct the synthetic profiles shown in Figure 4. The x-axis has units of logarithmic centimeters. 
The symbols and best fit (solid) curve are the same as those shown in the bottom panel of Figure 5. 
The dashed curves correspond to the ±2a errors (±4.2 x 10^^ cm~^) on the O I column density. 
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Fig. 7. — These four panels show the derivation of the final value for the H I column density in the 
direction of PG 1259+593. For all spectra shown in the four panels, Gaussian fits are superimposed 
on the data points, a) WSRT-only data, b) WSRT data "observed" by a 9'.7 beam. The WSRT 
channel maps were multiplied by a Gaussian with FWHM = 9'.7, integrated, and converted back 
to brightness temperature, c) Single-dish spectrum taken with the Effelsberg 100-meter telescope. 
This clearly shows the three components: Complex C at —128 kms~^, the IV- Arch at ~ —54 
kms~^, and the Milky Way ISM at ~ —5 kms""*^. d) Final combined spectrum. 
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Fig. 8. — See caption on next page. 
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Fig. 8. — Adopted model results for the H I Ly(5, Lye, and Ly9 lines, assuming (D/H)compiex c = 
0.0. The FUSE SiC2 data are shown as thin histogrammed lines. The O I absorption in the vicinity 
of the H I lines is represented by the smooth solid lines. The total (Ol + Hl + Dl) absorption 
model is shown as thick dashed curves. Each panel spans a wavelength range corresponding to a 
velocity range of ±500 kms~^ centered on the H I rest velocity, as indicated by the velocity scales 
above each panel. The locations of the H I and D I components are denoted by the vertical tick 
marks above each spectrum. Terrestrial H I airglow occurs between -100 and +100 kms~^^ in each 
panel and decreases in strength with increasingly higher lying lines in the Lyman series. The model 
reproduces the observed absorption well, except at the velocities of the Complex C D I components 
(#1 and #2), for which N(D I) was set to zero in this figure. (Note the extra absorption in the 
blue wings of the observed absorption profiles.) In all other components, D/H has been set to 
1.5 X 10~^. The source of the non-zero flux near 937.4 A in the Lye profile is unknown. It does 
not appear in the SiCl data. The steep absorption walls at negative velocities in the higher order 
H I Lyman-series transitions constrain the b-values of the Complex CHI components (see text for 
details) . 
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Fig. 9. — A breakdown of the component structure in the Ly^ hne. The model results are shown as 
dashed lines. For Complex C, D/H has been set to zero. The smooth solid line indicates the fit to 
the O I absorption along the sight line. Top: H I absorption produced by Complex C (components 
1 and 2). Middle: H I and D I absorption produced by the IV Arch (components 3 and 4). Bottom: 
H I and D I absorption produced by the Galactic ISM and positive velocity IVC (components 5 — 8). 
Note that the IV Arch and Galactic ISM components contribute no absorption at the velocity of 
the D I lines in Complex C. 
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Fig. 10.— Model fits (dashed lines) to the H I and D I Lya (STIS) and Ly/3 (FUSE) absorption hnes. 
The Ly/3 profile shown is night-only data obtained with the FUSE LiFl channel. The velocities 
of the components in the model are indicated above the Lyct spectrum. The H I absorption is so 
strong that radiation damping wings are visible in both lines; the D I absorption is overwhelmed by 
the H I absorption. These lines were used as a consistency check on the validity of the adopted H I 
model for the sight line. The absorption lines marked "IGM" are due to intervening intergalactic 
clouds along the sight line (see Richter et al. 2003). The emission feature near +500 kms~^ in 
the Ly/3 spectrum is residual O I airglow in these night-only data; the emission from these lines is 
much stronger than the emission in the O I lines used in our study of the O I absorption along the 
sight line. The widths of the airglow lines in the FUSE spectrum, Av w 100 kms~^, are broader 
than in the STIS spectrum because of the larger aperture used for the FUSE observation. 
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Fig. 11. — A comparison of the observed H I 21 cm emission in the direction of PG 1259+593 
(thin histogrammed hne) and the adopted H I model described in the text (heavy smooth hne) . At 
velocities —150 < vlsr < —100 kms^^ (i.e., Complex C velocities), the data are the interferometric 
measurements described in §2.3. At velocities outside this range, the data are singlc-dish (9^7 
beam) Effelsbcrg measurements. The model has been scaled to brightness temperature with the 
assumption that the emission is optically thin. 
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Fig. 12. — See caption on next page. 
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Fig. 12. — Adopted model results for the H I, D I, and O I absorption along the sight line assuming 
(D/H)compicx c = 2.6 X 10~^. The absorption features produced by the Lyman-series lines of Ly6- 
LyK in the FUSE SiC2 channel are shown as histogrammed lines. The total (O I + H I + D I) 
absorption model is shown as thick dashed curves. The smooth solid line isolates the contribution 
to the fit produced by O I absorption along the sight line. The velocity scale applies to the H I lines. 
The velocities of the SHI and D I components discussed in the text are indicated above the LjS 
and LyK lines. The downward pointing arrow in the top panel indicates the position of residual 
absorption that is not reproduced by the model (see §7.1). The horizontal error bar above the 
Lye absorption in the second panel indicates the absorption region used to calculate the equivalent 
widths in the D 1 curve-of-growth analyses discussed in the text. 
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Fig. 13. — Same as Figure 12, except for data from the FUSE SiCl channel. 
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Fig. 14. — See caption on next page. 
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Fig. 14. — Adopted model results for the H I, D I, and O I absorption along the sight line after 
incorporation of a weak H I absorption feature at -190 kms~^. This model has (D/H)(;;;ompicx c = 
1.8 X 10~^. The absorption features produced by the Lyman-series lines of hyd -LjK in the FUSE 
SiC2 channel are shown as histogrammed lines. The total (O 1 + H 1 + D I) absorption model is 
shown as thick dashed curves. The smooth solid line isolates the contribution to the fit produced 
by O I absorption along the sight line. The velocity scale applies to the H I lines. The velocities 
of the SHI and D I components discussed in the text are indicated above the hyS and Ly/t lines. 
The downward pointing arrow in the top panel indicates the position of the residual absorption 
modeled by the inclusion of a weak high-velocity H 1 feature. 
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Fig. 15. — Same as Figure 14, except for data from the FUSE SiCl channel. 
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Fig. 16. — Curve of growth results for the D I absorption in Complex C between ulsr = —160 
kms~^ and wlsr = —110 kms~^. The x-axis has units of logarithmic inverse centimeters. Filled 
circles are measurements based on data from the SiC2 channel. Open circles are measurements 
based on data from the SiCl channel. The points for the two channels have been offset slightly 
(plus or minus 0.02 dex) from their nominal values in the horizontal direction for clarity. Error 
bars are \a estimates; limits are 3cr estimates. Top: Single-component COG fit derived from the 
equivalent widths listed in Table 8. Curves of growth with b- values of 10 and 20 kms~^ and logN 
= 15.22 are shown as dashed lines. Bottom: Double-component COG model based on the fitting 
parameters in Table 7 used to construct the profiles shown in Figures 12 and 13. The parameters 
for the column densities, widths, and velocity separation of the two components arc listed in the 
top left corner. The model reproduces the observed equivalent widths as well as the COG fit shown 
in the top panel. The dashed lines indicate the curves corresponding to ±2 kms"^ changes in the 
b-value of the stronger component. The shape of the COG is insensitive to the choice of b-value 
for the weaker component. 
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Fig. 17. — Values of D/H as a function of metallicity as measured by [Si/H] (open symbols) 
or [0/H] (filled symbols). The data for this plot are contained in Table 10. The symbols are 
coded as follows: Complex C (triangle), high-redshift damped Lja systems (circles), high-redshift 
Lyman- limit systems (squares), local ISM (diamonds), and Galactic disk ISM (stars). Redshifts are 
indicated underneath the name of each extragalactic system. Error bars are la estimates from the 
original sources. The region between the horizontal dashed lines indicates the range of primordial 
D/H values consistent with recent cosmic microwave background measurements made by WMAP. 
The solid curve indicates the expected behavior of D/H as a function of metallicity for a simple 
chemical evolution model in which (D/H)p = 2.6 x 10^^ (see text). 



-66- 



Tablc 1. FUSE Observations of PG 1259+593^ 



Dataset 


Date 


# Exposures 


LiFl / SiCl 


LiF2 / SiC2 


Aperture 


Note^ 




(U.T. Start) 




ttot 

(ks) 


tngt 

(ks) 


ttot 

(ks) 


tngt 

(ks) 






P1080101 


2000-Feb-25 


12 


52.4 


14.2 


52.4 


14.2 


30" X 30" 




P1080102 


2000-Dec-25 


21 


56.8 


0.0 


0.0 


0.0 


30" X 30" 


1 


P1080103 


2001-Jan-29 


31 


52.0 


35.4 


82.1 


56.5 


30" X 30" 


2 


P1080104 


2001-Mar-12 


35 


105.9 


70.0 


103.2 


69.5 


30" X 30" 




P1080105 


2001-Mar-14 


35 


104.0 


70.1 


95.9 


69.8 


30" X 30" 




P1080106 


2001-Mar-17 


35 


63.6 


43.8 


66.7 


45.9 


30" X 30" 


3 


P1080107 


2001-Mar-19 


38 


95.7 


70.2 


96.1 


70.2 


30" X 30" 




P1080108 


2001-Mar-22 


13 


33.5 


24.1 


33.5 


24.1 


30" X 30" 




P1080109 


2001-Mar-28 


16 


32.7 


25.1 


32.0 


25.1 


30" X 30" 





'^Entries in this table include the dataset identification, U.T. date at the start of the 
observation, exposure times (total and night-only) in kiloseconds for the LiFl/SiCl channels 
and LiF2/SiC2 channels, and apertures used (30" x 30" = LWRS apertures). Exposure times 
are totals after screening for valid data with event bursts removed. The exposure times have 
been reduced when necessary to account for situations where the detector high voltage was 
down. These exposure times do not necessarily reflect the exact amount of time that the 
light of PG 1259+593 was in the aperture, except for the LiFl channel used for guiding. 

^Notes: 

1) SiCl, SiC2, and LiF2 were not aligned with LiFl, so no useful data were obtained in 

these channels. 

2) Detector 1 high voltage was down during exposures 18-29. 

3) High voltage was down on both detectors during exposures 6-17. High voltage was also 
down during exposure 22 on detector 1. 
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Table 2. HST/STIS Observations of PG 1259+593^ 



Dataset 


Date 


ttot 


Grating 


Aperture 




(U.T. Start) 


(ks) 






O63G05010-O63G05060 


2001-Jan-17 


14.4 


E140M 


0.2" X 0.06" 


O63G06010-O63G06060 


2001-Jan-17 


14.4 


E140M 


0.2" X 0.06" 


O63G07010-O63G07060 


2001-Jan-18 


14.4 


E140M 


0.2" X 0.06" 


O63G08010-O63G08060 


2001-Jan-18 


14.4 


E140M 


0.2" X 0.06" 


O63G09010-O63G09060 


2001-Jan-19 


14.4 


E140M 


0.2" X 0.06" 


O63G11010-O63G11040 


2001-Jan-19 


9.1 


E140M 


0.2" X 0.06" 


O63G10010-O63G10060 


2001-Dec-19 


14.4 


E140M 


0.2" X 0.06" 



^Entries in this table include the dataset identifications, U.T. date at 
the start of the observation, total exposure time for the datasets listed, 
grating, and aperture used for the observation. 
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Table 3. O I Lines' 



Wavelength 


log[/A(A)] 


Instrument 


S/Nb 


Used in 


Note'^ 


(A) 






(At; = 20 kms"^) 


I Fit? 




1302.168 


1.796 


STIS/E140M 


21 


Yes 


1 


1039.230 


0.974 


FUbE/LiF 


25 / 18 


Yes 




976.448 


0.509 


FUSE/biC 


13 / 11 


Yes 


2 


974.070 


-1.817 


FUSE/SiC 


13 / 11 


No 


3 


972.143 


-0.511 


FUSE/SiC 


13 / 11 


No 


4 


972.142 


-1.521 


FUbE/biC 


13 / 11 


No 


4 


971.738 


1.052 


FUSE/SiC 


13 / 11 


No 


5 


971.738 


0.304 


FUSE/biC 


13 / 11 


No 


5 


971.737 


0.872 


FUbE/biC 


13 / 11 


No 


5 


950.885 


0.176 


FUbE/biC 


13 / 10 


Yes 




948.686 


0.778 


FUbE/SiC 


13 / 10 


Yes 




937.840 


-0.085 


FUSE/SiC 


13 / 10 


No 


6 


936.630 


0.534 


FUSE/SiC 


13 / 10 


Yes 




930.257 


0.301 


FUSE/SiC 


13 / 11 


No 


7 


929.517 


0.329 


FUSE/SiC 


13 / 11 


Yes 




925.446 


-0.484 


FUSE/SiC 


13 / 11 


Yes 


8 


924.950 


0.155 


FUSE/SiC 


13 / 11 


Yes 




922.200 


-0.645 


FUSE/SiC 


13 / 11 


Yes 




921.857 


-0.001 


FUSE/SiC 


13 / 11 


Yes 




919.917 


-0.788 


FUSE/SiC 


13 / 11 


No 


9 


919.658 


-0.137 


FUSE/SiC 


13 / 11 


No 


10 



''This table lists prominent O 1 hnes in the FUSE and HST/STIS bandpasses. Lines 
with strengths log /A < —2.0 have been omitted from the table since they produce 
negligible absorption features. Wavelengths and /-values are from the atomic data 
compilations by Morton (1991, 2003). 

''Signal-to- noise ratio per 20 kms^^ bin in the continuum near the line. For the 
FUSE data, values are listed for either the LiFl / LiF2 channels (A > 1000 A) or the 
SiCl / SiC2 channels (A < 1000 A). 

'^Notes: 

(1) High-resolution spectrum; P II A1301.874 at 68 kms"-'^ with respect to O I 
features; O I terrestrial airglow emission at vlsr ~ 20.5 kms~^. 

(2) Partial blend at low velocities with high velocity C III A977.020. 

(3) Negligible absorption produced, blend with IGM line. 

(4) Blend with H 1 and D 1 Lyy. The published O I /-values for the resonance lines 
in this multiplet (multiplet #9) are probably a factor of at least 5 times too high (see 
Hebrard & Moos 2003). We do not use these lines in our analysis. 

(5) Possible partial blend with D I hyy. 

(6) Blend with H 1 Lye. 

(7) Blend with H 1 and D 1 LyC. 

(8) Close to, but does not blend with, D I Lyq. 

(9) Partial blend with H I LyK. 
(10) Blend with H I LyK. 
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Table 4. O I Velocity Model^ 



Component 


VLSR 
(kms~^) 


b 

(km s""-^) 


log [N (cm-2)] 


Note'' 


Complex C 


1 
2 

Total [1 2] 


-129.5 
-112.5 


6.0 + 1.0 
9.8+1;° 


15.85 + 0.15 
14.30 + 0.04 

15.8(3 ± 0.15 


1 
1 

2 


IV Arch 


3 
4 
5 

Total [3-5] 


-81.9 
-54.5 
-29.0 


~8 
-10 
-10 


15.08 
16.00 
14.78 
16.07 + 0.10 


3 
1 
3 
2 


Low Velocity ISM 


6 
7 

Total [6-7] 


-2.5 
+21.0 


—8 
-12 


16.08 
15.00 
16.11 + 0.10 


1 
1 
2 


Positive IVC 


8 


+68.9 


7.7 + 2.8 


13.43 + 0.08 


4 



^Uncertainties on all quantities in this table are la estimates. 
''Notes: 

1) Obvious component in high-resolution HST/STIS E140M 
metal-line data. 

2) Obvious component group in H I 21 cm emission. 

3) Weak component required to improve fit to the O I lines. 

4) Weak component visible in the O I A1302.168 line. Small 
inflection in O I A1039.230 line also present. 
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Table 5. Equivalent Widths of O I Lines in Complex C 



A 


log [/A(A)] 


This Paper 


Previous Work 






l^A(SiCl) 


VFA(SiC2) 


Richter ct al. (2001b) 


Collins et al. (2003) 


(A) 




(mA) 


(mA) 


(mA) 


(mA) 


1302.168 


1.796 


178 ± lO'' 


214 ± 20^ 


203 ± 8^ 


1039.230 


0.974 


90 ±18 


104 ± 12 


100 ± 13 


97 ±6 


948.685 


0.778 


60± 13 


82± 13 


84 ±25 


71 ± 12 


936.630 


0.534 


69 ±17 


66 ±10 


70 ±11 


73 ±10 


929.517 


0.329 


42± 17 


43± 11 


< 48 


54± 10 


924.950 


0.155 


40± 17 


36 ±10 




55 ±11 


921.857 


-0.001 


45 ±17 


38 ±13 







'^Errors are la uncertainties. Limits arc 3a estimates. The integration range in all cases was 
— 160 < I'LSR ^ —110 kms~^, which includes both Complex C components. 

'^Equivalent width measurement is from HST/STIS E140M data. 
^Equivalent width measurement is from partial set of HST/STIS E140M data. 
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Table 6. H I Velocity Model' 



Component 




b 


log [N (cm-2)] 


Note*^ 




(kms~^) 


(kms"-*^) 






Complex C 


1 


-127.4 


9.6 + 1.2 


19.94 + 0.06 


1 


2 


-112.5 


24 + 2 


18.43 + 0.30 


2 


Total [1-2] 






19.95 + 0.06 


3 


IV Arch 


3 


-81.9 




18.34 


4 


4 


-53.5 


~10 


19.43 


4 


5 


-29.0 


~10 


18.04 


4 


Total [3-5] 






19.48 + 0.01 


5 


Low Velocity ISM 


6 


-6.5 


~14 


19.63 


4 


7 


+21.0 


~12 


18.56 


4 


Total [6-7] 






19.67 + 0.02 


5 


Positive JVC 


8 


+68.9 


13 + 2 


16.69 + 0.08 


6, 7 



^Uncertainties on all quantities in this table are la estimates. 
^Notes: 

1) Component width and velocity set by interferometric 21 cm 
emission data. 

2) Component width constrained by the negative velocity wings 
of the H I absorption profiles. The large width implies that this 
component is likely to be a blend of several narrower unresolved 
sub-components. 

3) Total H I column density of Complex C set by interferomet- 
ric 21 cm data. 

4) Component width and velocity similar to the values for O I. 
Values constrained by 21 cm data. 

5) Total H I column density of IV Arch and ISM set by single- 
dish {9'.7 beam) Effelsberg 21 cm emission data. 

6) Component width set by the positive velocity wings of 
the H I absorption profiles. Component velocity set by O I 
A1302.168 line. 

7) Column density of positive IVC set by value of N(0 I) as- 
suming (O/H) = (0/11)0; no useful constraint on N(II I) is set 
by the H I Lyman-series lines. 
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Table 7. D I Velocity ModeP'*' 



Component 


•yLSR b lo^ 


r [N (cm-2)] 


Note'^ 




(kms~^) (kms~^) 








Complex C 






(without additional residual H I at 


-190 kms"^) 




1 


-127.4 8 ± 1 


15.36 


1 


2 


-112.5 15 


13.85 


2 


Total [1-2] 


15.37 ±0.10 


3 




Complex C 






(with additional residual H I at - 


190 kms~^) 




1 


-127.4 8 ± 1 


15.20 


1 


2 


-112.5 15 


13.69 


2 


Total [1-2] 


15.21 ±0.12 


4 


IV Arch 


3 


-81.9 8 


13.52 


5 


4 


-53.5 10 


14.61 


5 


5 


-29.0 10 


13.22 


5 


Total [3-5] 




14.66 


5 


Low Velocity ISM 


6 


-6.5 14 


14.81 


5 


7 


+21.0 12 


13.74 


5 


Total [6-7] 




14.85 


5 


Positive IVC 


8 


+68.9 13 


11.87 


5 



'^Velocities and widths of the D I absorption components 
were set equal to the values for H I in Table 6 unless noted 
otherwise. 

^Two sets of parameters are given for the Complex C ab- 
sorption. The first set is the fit without additional, weak high- 
velocity H I that could account for the discrepancy in the fits 
to the Lj5 and Lye lines shown in Figures 12 and 13. The 
second set includes a weak H 1 feature with N(II I) = 3 x 10^^ 
cm~^ and bm = 10 kms~^ in the fit, as shown in Figures 14 
and 15. 

'^Notes appear on next page. 
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^Notes: 

1) Component width allowed to vary between boi = 6 kms~^ and bni = 10 kms"^. 

2) Component width poorly constrained by data, but does not affect total column density of D I 
in Complex C. 

3) D/H ratio allowed to vary such that (D/H)compi = (D/H)comp2- Best fit is (D/H) = 2.6 x lO^^. 

4) D/H ratio allowed to vary such that (D/H)compi = (D/H)comp2- Best fit is (D/H) = 1.8 x 10~^. 

5) D/H ratio set equal to 1.5 x 10~^. 
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Table 8. Equivalent Widths of D I Lines in Complex C 



Line^ 


A 


log [/A(l)] 


WA(SiCl)^ 


WA(SiC2)b 




(A) 




(mA) 


(mA) 


LyS 


949.485 


1.122 


80 ±13 


111 ± 11 


Lye 


937.548 


0.865 


52± 10 


59 ±09 


Ly?? 


925.974 


0.470 


39 ±10 


42 ±09 


Ly6 


922.898 


0.311 


< 39 


< 32 


Lyi 


920.711 


0.170 


< 36 


< 35 



'^The Ly(^ line is omitted from this table because of 
substantial blending with O I absorption (see Figures 12 
and 13). 

''Errors are la uncertainties. Limits are 3cr estimates. 
The integration range in all cases was —160 < vlsr < 
— 110 kms~^, which includes both Complex C compo- 
nents. 
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Table 9. Summary of Complex C Results 



Quantity 




Value lb lo" 


Value lb 2a 






(68% Conf.) 


(95% Conf.) 


N(H I) (cm- 




(9.0 lb 1.0) X 10^9 


(9.0 lb 1.6) X 10^9 


N(D I) (cm- 




(2.0 zb 0.6) X 10^^ 


(2.0 lb 0.9) X 10^5 


N(0 I) (cm- 




(7.2 lb 2.1) X 10^5 


(7.2 lb 4.2) X 10^5 


D/H 




(2.2 lb 0.7) X 10"^ 


(2.2 lb 1.1) X 10-^ 


0/H 




(8.0 lb 2.5) X 10^5 


(8.0 lb 4.2) X 10"^ 


D/0 




0.28 lb 0.12 


0.28 lb 0.20 



Note - Errors and confidence limits include statistical 
uncertainties as well as systematic uncertainties discussed 
in the text. 
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Table 10. Light Element Abundance Ratios in Different Environments^ 



System 


z or 
[d (pc)] 


log N(HI [cm-2]) 


[O/H]" 


D/H 
(10-5) 


D/0 


Reference 


High Velocity Gas 


Complex C 


[> 5000] 


19.95 


-0.79+8:}^ 


2.2 + 0.7 


0.28 + 0.12 


This paper 


Milky Way ISM 


LISM (7 l.o.s.) 


[< 180] 


17.93 - 20.14 


-o.2i±8:8i 


1.52 + 0.08 


0.040 + 0.002 


Moos et al. (2002)"= 


7^ Vel 


[258] 


19.71 


(-0.23) 


2.i8+8:?i 




Sonneborn et al. (2000)=^' 


CPup 


[430] 


19.96 




i.42+8:it 




Sonneborn et al. (2000)'' 


<5 Ori A 


[500] 


20.19 


-o.24±8:g? 


o.74+8:Ji 


o.o26+8:88i 


Jenkins et al. (1999)"-'' 


HD 195965 


[794] 


20.95 


+0.13+8:8! 


0.85+8:11 


o.oi3+8:88i 


Hoopes et al. (2003) 


HD 191877 


[2200] 


21.05 


-0.20+8:?! 


o.78+8:?i 


o.o25+8:8m 


Hoopes et al. (2003) 


Damped Lya Systems 


HS 0105+1619 


2.536 


19.42 


-1.73 + 0.02 


2.54 + 0.23 


2.8 + 0.3 


O'Meara et al. (2001) 


Q 0347-3819 


3.025 


20.63 


-0.68 + 0.06 


3.75 + 0.25 


0.37 + 0.03 


Levshakov et al. (2002) 


Q 2206-199 


2.076 


20.44 


(-2.23) 


1.65 + 0.35 




Pettini & Bowen (2001) 


Q 1243+3407 


2.526 


19.73 


-2.79 + 0.05 


2.42+8:ii 


30 + 3 


Kirkman et al. (2003) 


Lyman Limit Systems 


PKS 1937-1009 


3.572 


17.86 


(-2.26) 


3.25 + 0.30 




O'Moara et al. (2001) 


Q 0130-403 


2.799 


16.66 


( 2.6) 


< 6.8 




O'Meara et al. (2001) 


Q 1009+299 


2.504 


17.39 


(-2.53) 


3.98+8:1? 




O'Meara et al. (2001) 



^Uncertainties arc Icr estimates. 

''[O/H] = log N(0 1)/N(H 1) - log (0/H)©, where (O/H)© = 4.90 x 10"* (AUende Prieto et al. 2001). All values of [O/H] 
have been scaled to this reference abundance. Items in parentheses are estimates of [Si/H] when no estimate of [O/H] is 
available. These values of [Si/H] should only be considered approximate substitutions for [O/H] since ionization corrections 

for Si may bo very large, especially for the Lyman-limit systems. 

'^Values are [O/H], D/H, and D/O are weighted averages for the seven sight lines considered. Measurements for each sight 
line can be found in Kruk et al. (2002), Friedman et al. (2002), Sonneborn et al. (2002), Lemoine et al. (2002), Lehner et 
al. (2002), Wood ct al. (2002), and Hcbrard ct al. (2002). Additional measurements of D/O can be found in Hebrard & 
Moos (2003), who measure a local ISM value of D/O = 0.038 + 0.002. 

'*N(Si II) is from Fitzpatrick & Spitzer (1994). 

"Errors on D/H have been converted from 90% confidence to lu estimates (see Moos et al. 2002). 
'N(0 1) is from Meyer et al. (1998). 
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Table 11. Searching for Deuterium in Other Complex C Sight Lines 



Sight Line 


I 

(°) 


b 

(°) 


Q3. 

(°) 


Comment 


PG 1626+554 


84.51 


42.19 


27.7 


Velocity structure similar to PG 1259+593, but flux is lower and 
Complex C component is less pronounced than Milky Way ISM. 


3C351 


90.08 


36.38 


29.5 


Ly-limit system at z ^ 0.22 attenuates quasar flux shortward 
of ^ 1115 A. 


Mrk 290 


91.49 


47.95 


20.0 


Potentially good velocity structure, but quasar is faint. 


Mrk 876 


98.27 


40.38 


22.7 


Good Complex C structure, but D I absorption blends with H I. 


Mrk 817 


100.30 


53.48 


12.2 


Complex C at relatively low velocity (-109 kms~"^), 
leading to blending of D I and H I. 


PG 1351+640 


111.89 


52.02 


7.8 


Strong intermediate velocity structure confuses Complex C 
D I absorption. 


Mrk 279 


115.04 


56.86 


3.2 


Complicated structure at intermediate velocities blends with 
Complex C absorption. 



Angular distance from PG 1259+593 sight line. 



